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Abstract - We discuss the consequences of the introduction of a quantum of 
time To in the formalism of non-relativistic quantum mechanics, by referring ourselves, in 
^ ! particular, to the theory of the chronon as proposed by P.Caldirola. Such an interesting 

"finite difference" theory, forwards — at the classical level — a solution for the motion 
of a particle endowed with a non-negligible charge in an external electromagnètic field, 
overcoming all the known difficulties met by Abraham-Lorentz's and Dirac's approaches 
(and even allowing a clear answer to the question whether a free falling charged particle 
qh| does or does not emit radiation), and — at the quantum level — yields a remarkable mass 

spectrum for leptons. 

After having briefly reviewed Caldirola's approach, our first aim is to work out, 
discuss, and compare to one another the new representations of Quantum Mechanics 
(QM) resulting from it, in the Schròdinger, Heisenberg and density-operator (Liouville- 
von Neumann) pictures, respectively. 

Moreover, for each representation, three (retarded, symmetric and advanced) for- 
mulations are possible, which refer either to times t and t — tq, or to times t — to/2 
and t + To/2, or to times t and t + tq, respectively. It is interesting to notice that, 
when the chronon tends to zero, the ordinary QM is obtained as the limiting case of 
the "symmetric" formulation only; while the "retarded" one does naturally appear to 
describe QM with friction, i.e., to describe dissipative quantum systems (like a particle 
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moving in an absorbing medium). In this sense, discretized QM is much richer than the 
ordinary one. 

We also obtain the (retarded) finite-difference Schròdinger equation within the Feyn- 
man path integral approach, and study some of its relevant solutions. We then derive 
the time-evolution operators of this discrete theory, and use them to get the finite- 
difference Heisenberg equations. 

When discussing the mutual compatibility of the various pictures listed above, we 
find that they can be written down in a form such that they result to be equivalent (as 
it happens in the "continuous" case of ordinary QM), even if the Heisenberg picture 
cannot be derived by "discretizing" directly the ordinary Heisenberg representation. 

Afterwards, some typical applications and examples are studied, as the free particle, 
the harmònic oscillator and the hydrogen atom; and various cases are pointed out, for 
which the predictions of discrete QM differ from those expected from "continuous" QM. 

At last, the density matrix formalism is applied to the solution of the measurement 
problem in QM, with very interesting results, as for instance a natural explication of 
"decoherence" , which reveal the power of dicretized (in particular, retarded) QM. 
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1 Introduction 



The idea of a discrete temporal evolution is not a new one and, as almost all the 
physical ideas, has from time to time been recovered from oblivionjj For instance, in 
classical Greece this idea came to light as part of the atomistic thought. In the Middle 
Age, belief in the discontinuous character of time was at the basis of the "theistic 
atomism" held by the Arabic thinkers of the Kalamj2]. In Europe, discussions about 
the discreteness of space and time can be found for example in the writings of Isidore of 
Sevilla, Nicolaus Boneti and Henry of Harclay, who discussed the nature of continuum. 
In more recent times, the idea of the existence of a fundamental interval of time was 
rejected by Leibniz, since it was incompatible with his rationalistic philosophy. Within 
modern physics, however, Planck's famous work on black body radiation inspired a new 
view of the subject. In fact, the introduction of the quanta opened a wide range of new 
scientific possibilities regarding the way the physical world can be conceived. Including 
considerations, like those in the present paper, on the "discretization" of time within 
the framework of Quantum Mechanics (QM). 

In the early years of our century, Mach regarded the concept of continuum to be a 
consequence of our physiological limitations: <<...le temps et l'espace ne représentent, 
au point de vue physiologique, qu'un continue apparent, qu'ils se composent très 
vraisemblablement d'elements discontinus, mais qu'on ne peut distinguer nettement 
les uns des autres.>>[3] Also Poincaré took into consideration the possible existence 
of what he called an "atom of time": the minimum amount of time which allows to 
distinguish between two states of a system. [1] Finally, in the twenties, J.J.Thomson[5j 
suggested the elèctric force to act in a discontinuous way, producing finite increments 
of momentum separated by finite intervals of time. Such a seminal work has ever since 
inspired a series of papers on the existence of a fundamental interval of time, named 
chronon; even if the repercussions of all that work was small, at that time. A further 
seminal article was the one by Ambarzumian and Ivanenko[6], appeared in 1930, which 
assumed space-time as being discrete and also stimulated a large number of subsequent 
papers. 

It is important to stress that, in principle, time discretization can be introduced in 
two distinct (and completely different) ways: 



*Historical aspects related to the introduction of a fundamental interval of time in Physics can be 
found in F.Casagrande, ref.pQ. 
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(1) by attributing to time a discrete structure, i.e., by regarding time not as a 
continuum, but as a one-dimensional "lattice" ; 



(2) by considering time as a continuum, in which events can take place (discontinu- 
ously) only at discrete instants of time. 

Almost all the attempts to introduce a discretization of time followed the flrst way, 
generally as part of a more extended procedure in which the space-time as a whole is 
considered as intrinsically discrete (four-dimensional lattice). Recently, also T.D.Lee[7] 
introduced a time discretization on the basis of the finite number of experimental mea- 
surements performable in any finite interval of time. For an early approach in this 
direction, see e.g. T.Tati[8j and references therein. 

The second approach was first adopted in the twenties (e.g., by LeviJTÜ] and by 
Pokrowski|7Uj). after Thomson's work, and resulted in the first real example of a 
theory based on the existence of a fundamental interval of time: the one set forth by 
Caldirola, in the fifties[TÜl [TT]. Namely, Caldirola formulated a theory for the classical 
electron, with the aim of providing a consistent (classical) theory for its motion in 
an electromagnètic field. In the late seventies, Caldirola extended its procedure to 
non-relativistic quantum mechanics. 

It is known that the classical theory of the electron in an electromagnètic field 
(despite the efforts by Abraham,[T2] Lorentz.[13| Poincaré[Il] and Dirac,[T5] among 
others) actually presents many serious problems; except — of course — when the 



finíte-dífference equations, Caldirola developed a theory in which the main difficulties 
of Dirac's theory were overcome. As we shall see, in his relativistically invariant 
formalism the chronon characterizes the changes suffered by the dynamical state of the 
electron when it is submitted to external forces. So that the electron will be regarded 
as an (extended-like) object, which is point-like only at discrete positions x n (along 
its trajectory) such that the electron takes a quantum of proper time to travel from 
one position to the following one [or, rather, two chronons: see the following]. It is 

^It is interesting to note that all those problems have been — necessarily — tackled by Yaghjian in 
his recent book.[T6] when he faced the question of the relativistic motion of a charged, macroscopic 
sphere in an external electromagnètic field. 




By replacing Dirac's 



differential equation by two 
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tempting to examine extensively the generalization of such a theory to the quantum 
domain; and this will be performed in the present work. Let us recali that one of the 
most interesting aspects of the discretized Schròdinger equations is that the mass of 
the muon and of the tau lepton followed as corresponding to the two levels of the first 
(degenerate) excited state of the electron. 

In conventional QM there is a perfect equivalence among its various pictures: 
Schròdinger's, Heisenberg's, density matrix's . When discretizing the evolution equa- 
tions, we shall succeed in writing down those pictures in a form such that they result to 
be still equivalent. However, in order to be compatible with the Schròdinger representa- 
tion, our Heisenberg equations cannot, in general, be obtained by a direct discret ization 
of the continuous Heisenberg equation. 

The plan of this work is the following. In Chaptj2] we present a brief review of the 
main classical theories of the electron, including Caldirola's. In ChaptElwe introduce 
the three discretized f ornis ("retarded", "advanced", "symmetrical") of the Schròdinger 
equation, analyze the main characteristics of such formulations, and derive the retarded 
one from Feynman's path integral approach. In ChaptHJ our discrete theory is applied 
to some simple quantum systems, such as the harmònic oscillator, the free particle 
and the hydrogen atom. The possible experimental deviations from the predictions of 
ordinary QM are investigated. In Chaptl5l a new derivation of the discretized Liouville- 
von Neumann equation, starting from the coarse graíned hypothesis, is presented. Such 
a representation is then adopted for tackling the measurement problem in QM, with 
rather interesting results. Finally, a discussion on the possible interpretation of our 
discretized equations can be found in Chaptjül 

2 The Introduction of the Chronon in the Classical 
Theory of the Electron 

Almost a century after its discovery, the electron continues to be an object waiting 
for a convincing description (cf., e.g., [32]), both in classical and quantum electrody- 
namics. As Schròdinger put it, the electron is — still — a stranger in electrodynamics. 
Maxwell's electromagnetism is a field theoretical approach in which no reference is made 
to the existence of material corpuscles. Thus, one may say that one of the most con- 
troversial qüestions of the twentieth century physics, the wave-particle paradox, is not 
characteristic of QM only. In the electron classical theory, matching the description 



5 



of the electromagnètic fields (obeying Maxwell equations) with the existence of charge 
carriers like the electron is still a challenging task. 

The hypothesis that elèctric currents could be associated with charge carriers was 
already present in the early "particle electrodynamics" formulated in 1846 by Fechner 
and Weber. But such an idea was taken into consideration again only a few dècades 
later, in 1881, by Helmholtz. Up to that time, electrodynamics had been developed 
on the hypothesis of an electromagnètic continuum^E\ and of an ether. In that same 
year, J.J.Thomson wrote his seminal paper in which the electron mass was regarded as 
purely electromagnètic in nature. Namely, the energy and momentum associated with 
the (electromagnètic) fields produced by an electron were held entirely responsible for 
the energy and momentum of the electron itself . [19] 

Lorentz's electrodynamics, which described the particle-particle interaction via elec- 
tromagnètic fields by the famous force law 



p being the charge density of the particle on which the fields act, dates back to the 
beginning of the 1890 decade. The electron was finally discovered by Thomson in 1897, 
and in the following years various theories appeared. The famous (pre-relativistic) 
theories by Abraham, Lorentz and Poincaré regarded it as an extended-type object, 
endowed again with a purely electromagnètic mass. As well-known, in 1903 Abraham 
proposed the simple-minded (and qüestionable) model of a rigid sphere, with a uniform 
elèctric charge density on its surface. The theory of Lorentz (1904) was quite similar, 
trying to improve the situation with the mere introduction of the effects resulting from 
the Lorentz-Fitzgerald contraction. 

2.1 The theory of the electron by Abraham— Lorentz 

A main difficulty for an accurate description of the electron motion was the inclusion of 
the radiation reaction, i.e., of the effect produced on such a motion by the fields radiated 
by the particle itself. In the model proposed by Abraham-Lorentz the assumption of a 
purely electromagnètic structure for the electron implied that 



where F p is the self-force due to the self-fields of the particle, and F ex t is the external 
force. According to Lorentz's law, the self-force was given by 





F p + F oxt — , 



(2) 
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F 

r p 



J p (E p + iv A B p ) d 3 r 



where E p and B p are the fields produced by the charge density p itself, according to 
the Maxwell-Lorentz equations. For the radiation reaction force, Lorentz obtained the 
following expression: 

where = (47reo) _1 [in the following, whenever convenient, we shall assume units such 
that numerically k = 1], and where 

i //4£MïO d 3 rd 3 r , 



2 

is the electrostatic self-energy of the considered charge distribution, and R is the radius 
of the electron. All the terms in the sum are structure dependent; i.e., they depend 
on R and on the charge distribution. By identifying the electromagnètic mass of the 
particle with its electrostatic self-energy 



m el = — , 



it was possible to write eq.© as 

4 



3 m e iv -T = F ext , (4) 
so that one got: 

r=~à(i + o(fl)) , (5) 

which was the equation of motion in the Abraham-Lorentz model. Quantity T is the 
radiation reaction force, namely, the reaction force acting on the electron. A problem 
with equation (4) was constituted by the factor g. In fact, if the mass is supposed to 
be of electromagnètic origin only, then the total momentum of the electron would be 
given by 

P = 3— v , (6) 

which is not invariant under Lorentz transformations. That model, therefore, was non- 
relativistic. Finally, we can observe from equation ([3]) that the structure dependent 
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terms are functions of higher derivatives of the acceleration; even more, the resulting 
differential equation is of the thírd order, so that initial position and initial velocity are 
not enough to single out a solution. In order to suppress the structure terms, one ought 
to reduce the electron to a point, (R — > 0), but in this case the self-energy W e \ and 
mass m el would diverge! 

After the emergence of the special theory of relativity or, rather, after the publication 
by Lorentz in 1904 of his famous transformations, some attempts were made to adapt 
the model to the new requirements. [2Ü1 I2T1 122] Abraham himself (1905) succeeded in 
deriving the following gener alization of the radiation reaction term ([5]): 



A solution for the problem of the electron momentum non-covariance was proposed 
by Poincaré in 1905, by the addition of cohesive forces of non-electromagnetic character, 
which — however — made the electron no longer purely electromagnètic in nature. 

On the other hand, electrons could not be considered pointlike, due to the obvious 
divergence of their energy when R — > 0; thus, a description of the electron motion 
could not forget about the structure terms. Only Fermi succeeded to show, later, that 
the correct relation for the momentum of a purely electromagnètic electron could be 
obtained without Poincaré's cohesive forces. [23] 

2.2 Dirac's theory of the classical electron 

Notwithstanding its inconsistències, the theory by Abraham-Lorentz was the most ac- 
cepted theory of the electron, until the publication of Dirac's theory in 1938. During 
the long period in between, as well as afterwards, however, various further attempts to 
solve the problem were set forth, either by means of extended-type models (Mie, Page, 
Schott, etc.[24j), or by trying again to treat the electron as a pointlike particle (Fokker, 
Wentzel, etc.[25]). 

Dirac's approach[T5] is the best known attempt to describe the classical electron. 
It by-passed the critical problem of the previous theories of Abraham and Lorentz by 
working out for the pointlike electron a trick which avoided divergences! By using the 
conservation laws of energy and momentum, and Maxwell equations, Dirac calculated 
the flux of the energy-momentum 4-vector through a tube of radius e < iü (quantity 
R being the radius of the electron at rest) surrounding the world line of the particle, 
and obtained: 




(7) 
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m^ = F, + T l·ll (8) 

where is the Abraham 4- vector ((7|), that is, the reaction force acting on the electron 
itself; and is the 4-vector that represents the external field acting on the particle: 

F„ = -F^vï . (9) 

According to such a model, the rest mass rriQ of the electron is the limiting, finite 
value obtained as the difference of two quantities tending to infinity when R — > 0: 




m 

the procedure followed by Dirac being an early example of elimination of divergences 
by means of a subtractive method. 

At the non-relativistic limit, Dirac's equation goes into the one previously obtained by 
Abraham-Lorentz: 

dv 2e 2 d 2 v / 1 



ro °dF-3?dF = e l E + - vAB ) • (10) 



C 



except for the fact that in Abraham-Lorentz's approach diverged. The last equation 
shows that the reaction force equals | ^ . 

Dirac's dynamical equation (jSj) was later reobtained also from different, improved 
models. [26J Wheeler and Feynman, for example, rederived eq.© by basing electromag- 
netism on an action principle applied to partides only, via their absorber hypothesis. |27j 
However, eq. ([8"]) also presents many troubles, related to the infinite many non-physical 
solutions that it possesses. Actually, as already mentioned, it is a third-order differen- 
tial equation, requiring three initial conditions for singling out one of its solutions. In 
the description of a free electron, e.g., it even yields "self-accelerating" solutions (run- 
away solutions), [28] for which velocity and acceleration increase spontaneously and 
indefmitely. [25J Moreover, for an electron submitted to an electromagnètic pulse, fur- 
ther non-physical solutions appear, related this time to pre-accelerations^29\. If the 
electron comes from infinity with a uniform velocity vq and, at a certain instant of 
time t , is submitted to an electromagnètic pulse, then it starts accelerating before t . 
Drawbacks like these motivate further attempts to find out a coherent model for the 
classical electron. 
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2.3 Caldirola's theory for the classical electron 

Among the various attempts to formulate a more satisfactory theory, we want to focus 
our attention on the one proposed by P.Caldirola. Like Dirac's, Caldirola's theory is 
also Lorentz invariant. Continuity, in fact, is not an assumption required by Lorentz 
invariance. [30] The theory postulates the existence of a universal interval r of proper 
time, even if time flows continuously as in the ordinary theory. When an external force 
acts on the electron, however, the reaction of the particle to the applied force is not 
continuous: The value of the electron velocity is supposed to jump from u^r — r ) 
to u M (r) only at certain positions s n along its world line; these discrete positions being 
such that the electron takes a time r for travelling from one position s n _i to the next 
one s n . 

In this theory [TÜ] the electron, in principle, is still considered as pointlike, but the 
Dirac relativistic equations for the classical radiating electron are replaced: (i) by a 
corresponding finite-difference (retarded) equation in the velocity u^ij) 

m j (r) u v (r) \ 
— < Uf, (r) - u M (r - T ) H ~ 2 [u v (r) - u v (r - r )J > = 

= e -F, v {T)u v {r), (11) 

which reduces to the Dirac equation (jSJ) when r — > 0, but cannot be derived from it (in 
the sense that it cannot be obtained by a simple discretization of the time derivatives 
appearing in Dirac's original equation); and (ii) by a second equation, connecting this 
time the "discrete positions" x^(t) along the world line of the particle. In fact, the 
dynamical law above is by itself unable to specify univocally the variables m m (t) and 
x fx{ T ) which describe the motion of the particle. Caldirola named it the transmission 
law. 

•Efj, (nr ) - x M [(n - 1) r Q ] = y {u^ (nr ) - n M [(n - 1) r ]} , (12) 

which is vàlid inside each discrete interval tq, and describes the internal or microscopic 
motion of the electron. 

In these equations, w m (t) is the ordinary four- vector velocity satisfying the condition 

w /t (r)'U ít (r) = — c 2 for r = nr 
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where n — 0, 1, 2, ... and — 0, 1, 2, 3; F**" is the external (retarded) electromagnètic 
field tensor, and the quantity 

^ = e = ~ 6.266 x 10~ 24 s (13) 

2 3 m c ó 

is defined as the chronon associated with the electron (as it will be justified below). 
The chronon 6q = Tq/2 depends on the particle (internal) properties: namely, on its 
charge e and rest mass m . 

As a result, the electron happens to appear eventually as an extended-like[3T] parti- 
cle, with internal structure, rather than as a pointlike object (as initially assumed). For 
instance, one may imagine that the particle does not react instant aneously to the action 
of an external force because of its finite extension (the numerical value of the chronon 
is of the same order as the time spent by light to travel along an electron classical 
diameter). As already said, eq. fllip describes the motion of an object that happens to 
be pointlike only at discrete positions s n along its trajectory[TÜ]; even if both position 
and velocity are still continuous and well-behaved functions of the parameter r, since 
they are differentiable functions of r. 

It is essential to notice that a discrete character is given to the electron merely by 
the introduction of the fundamental quantum of time, without any need of a "model" 
for the electron. Actually it is well-known that many difficulties are met not only by 
the strictly pointlike models, but also by the extended-type particle models ( "spheres" , 
"tops", "gyroscopes", etc). In A.O.Barut's words, "If a spinning particle is not quite a 
point particle, nor a sòlid three dimensional top, what can it be?" . We deem the answer 
stays with a third type of models, the "extended-like" ones, as the present theory; or as 
the (related) theoretical approach [3Tl [T7] in which the center of the pointlike charge is 
spatially distinct from the particle center-of-mass. Anyway, it is not necessary to recali 
that the worst troubles met in quantum field theory (e.g., in quantum electrodynamics), 
like the presence of divergències, are due to the pointlike character still attributed to 
(spinning) partides; since the problem of a suitable model for elementary partides was 
transported, unsolved, from classical to quantum physics. One could say that problem 
to be the most important in modern particle physics. 

Equations (ÍTTI) and (ÍT2|) together provide a full description of the motion of the 
electron. Notice that the global, "macroscopic" motion can be the same for different 
solutions of the transmission law. The behaviour of the electron under the action of 
external electromagnètic fields is completely described by its macroscopic motion. 
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As in Dirac's case, the equations above are invariant under Lorentz transformations; 
but, as we are going to see, are free from pre-accelerations, self-accelerating solutions, 
and problems (that had raised great debatés in the first half of the century) with the 
hyperbolic motion. 

In the non-relativistic limit the previous (retarded) equations reduced to the form 



which can be obtained — this time- from eq. (1101) . by directly replacing the time deriva- 
tives by the corresponding finite-difference expressions. The macroscopic equation (jT3J) 
had already been obtained also by other authors[2Tl [221 EU El] for the dynamics of 
extended-type electrons. 

The important point is that eqs. fíTTT) . fíT2|) . or eqs. fíT4l) . fíT5|) . allow to overcome the 
difficulties met with the Dirac classical equation (jSj). In fact, the electron macroscopic 
motion is completely determined once velocity and initial position are given. Solutions 
of the relativistic equations (ÍTT1) . (ÍT2|) for the radiating electron — or of the corresponding 
non-relativistic equations ÍÍT4"|) . ÍÍT51) — were obtained for several problems, the resulting 
motions never presenting unphysical behaviour; so that the following qüestions can be 
regarded as having been solved[TÜ]: A) exact relativistic solutions: 1) free electron 
motion; 2) electron under the action of an electromagnètic pulse;[35j 3) hyperbolic 
motion[32]; B) non-relativistic approximate solutions: 1) electron under the 
action of time-dependent forces; 2) electron in a constant, uniform magnètic field[36j; 
3) electron moving along a straight line under the action of an elàstic restoring for ce [37j. 

Before going on, let us briefly study the electron radiation properties as deduced 
from the finite-difference relativistic equations fi 111) . (I12p . with the aim of showing the 
advantages of the present formalism with respect to the Abraham-Lorentz-Dirac one. 
Such equations can be written[32| [TÜ] 



m 



[ v (t) - v (í - r )] = e E (t) + -v (t) A B (t) 



(14) 




(15) 



AQ M (r) 



+ ^(r) + ^(r) 



(16) 



where: 




(17) 
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Mr) = ~ \ UAT) p T) K(r + r ) + ^(r-r )-2^(r)] 



-W) = -^T \^—# [ M (r + r )-M (r-T )] \ . (19) 



In eq. (1161) . the first term AQq/tq represents the variation per unit of proper time (in 
the interval r — r to r) of the particle energy-momentum vector. The second one, 
-R M (r), is a dissipative term, since it contains only even derivatives of the velocity as 
one can prové by expanding u v {j + r ) and u v [j — r ) in terms of r ; furthermore, it 
is never negative, [TÜl [32] and can therefore represent the energy-momentum radiated 
by the electron in the unit of proper time. The third term, S^r), is conservative and 
represents the rate of change in proper time of the electron reaction energy-momentum. 
The time component (fi = 0) of eq. (ÍTü|) writes: 

T(r) - T(r - r ) + + = 

To 

quantity T(r) being the kinetic energy 

so that in eq. ( Í2DI) the first term replaces the proper-time derivative of the kinetic energy, 
the second one is the energy radiated by the electron in the unit of proper time, S (t) is 
the variation rate in proper time of the electron reaction energy (radiative correction), 
and P ext (r) is the work done by the external forces in the unit of proper time. 

We are now ready to show, e.g., that eq. fl2Ü|) yields a clear explanation for the 
origin of the so-called "acceleration energy" (Schott energy), appearing in the energy- 
conservation relation for the Dirac equation. In fact, by expanding in power series with 
respect to r the l.h. sides of eqs. ííTülKÍÍTTl) . (fTg]) . (fl9"]) for fj, — 0, and keeping only the 
first-order terms, one gets 

T(r) - T(r - r ) ^ dT _ 2 e*_ doç . 
To dr 3 c 2 dr ' 

Mr) * 7r=Fïïí aX; (23) 

Sq{t) ~ ; (24) 
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where a M is the four-acceleration 

drmu^ drmu^ 



= = 7 



drmr drmt 

quantity 7 being the Lorentz factor. Therefore eq. fl2Ü|) . to the first order in tq, becomes: 

or, passing from the proper time r to the observer's time t: 

dt 3 c 2 dí 3 c M v ; dí v ; 

The last relation is identical with the energy conservation law found by Fulton and 
Rohrlich|41j for the Dirac equation. In eq.(26) the derivative of (2e 2 /3c 2 )ao appears, 
which is just the acceleration energy. Our approach clearly shows that it arises only by 
expandíng in a power series of r the kinetic energy increment suffered by the electron 
during the fundamental proper-time interval To; whilst such a Schott energy (as well 
as higher-order energy terms) does not show up when adopting the full formalism of 
finite-difference equations. We shall come back to this important point in Subsection 
2.4. 



Let us finally observe[10j that, when setting 



""' [u,(r)u u (r - T ) - u,(r - t )u u (t)] = F^f , (27) 



ecT 

the relativistic equation of motion (ÍTT|) reads: 



e 
c 

confirming that F^} represents the (retarded) self-field associated with the moving 
electron. 



2.4 The three alternative formulations of Caldirola's theory 

Two more (alternative) formulations are possible of Caldirola's equations, based on 
different discretization procedures. In fact, equations (TTTT) and (fl2"|) describe an intrin- 
sically radiating particle. And, by expanding equation fíTTj) in terms of r , a radiation 
reaction term appears. Caldirola called those equations the retarded form of the electron 
equations of motion. 
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By rewriting the fmite-difference equations, on the contrary, in the form: 



— < (r + T ) - u l·l (r) + -£ — [«„ (r + r ) - u„ (r)J 



-F M1/ (r)^(r), (29) 

c 



.t m [(n + 1) r ] - x„ (nr ) = r n M (nr ) , (30) 

one gets the advanced formulation of the electron theory, since the motion — according to 
eqs. (I2"9~j) . (Í3"0l — is now determined by advanced actions. In contrast with the retarded 
formulation, the advanced one describes an electron which absorbs energy from the 
external world. 

Finally, by adding together retarded and advanced actions, Caldirola wrote down 
the symmetric formulation of the electron theory: 



J , , , f s . Up (r) u u (r) 
— <j (r + r ) - Uft [t - T ) + — — - 2 \u v (r + r ) - u v (r - r )J 



-F^(r)n,(r), (31) 

c 



[(n + 1) r Q ] - s M ((n - 1) r ) = 2t m m (nr ) , (32) 
which does not include any radiation reactions, and describes a non radiating electron. 

Before closing this brief introduction to Caldirola's theory, it is worthwhile to present 
two more relevant results derived from it, one of them following below, in the next 
Subsection. If we consider a free particle and look for the "infernal solutions" of the 
equation ( ÍTBl . we get — for a periodical solution of the type 

x = — p c sm 

y = —/3 c cos 
i = 



T 
Í2-KT 

V r 



(which describes a uniform circular motion) and by imposing the kinetic energy of the 
interna! rotational motion to equal the intrinsic energy moc 2 of the particle — that the 
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amplitude of the oscillations is given by ffi = §. Thus, the magnètic moment corre- 
sponding to this motion is exactly the anomalous magnètic moment of the electron, [38J 
obtained here in a purely classical context: 



1 e 



3 



Va. = -r- 



47r rriQC 2 

This shows that the anomalous magnètic moment is an intrinsically classical, and not 
quantum, result; and the absence of h in the last expression is a remar kable confirmation 
of this fact. 



2.5 Hyperbolic motions 

In a review paper on the theories of electron including radiation-reaction effects, 
Erber[3ü] criticized Caldirola's theory for its results in the case of hyperbolic motion. 

Let us recali that the opinion of Pauli and von Laue (among others) was that 
a charge performing uniformly accelerated motions — e.g., an electron in free fall — 
could not emit radiation. That opinion was strengthened by the invariance of Maxwell 
equations under the group of conformal transformations, [1Ü] which in particular includes 
transformations from rest to uniformly accelerated motions. Since the first dècades of 
the twentieth century, this had been — however — an open question, as the works by 
Born and Schott had on the contrary suggested a radiation emission in such a case. In 
1960, Fulton and Rohrlich^T] demonstrated that from Dirac's equation for the classical 
electron the emission of radiation during the hyperbolic motion follows. 

A solution of this paradox is possible within Caldirola's theory, and it was worked 
out in ref. [32] By analysing the energy conservation law for an electron submitted to an 
external force and following a procedure similar to that of Fulton and Rohrlich, Lanz 
obtained the expression (1201 above. By expanding it in terms of r and keeping only 
the first order terms, he arrived at expression (1251) . identical to the one obtained by 
Fulton and Rohrlich, in which — let us repeat — the Schott energy appears: a term 
that Fulton and Rohrlich (having obtained it from Dirac's expression for the radiation 
reaction) interpreted as a part of the internal energy of the charged particle. 

For the particular case of hyperbolic motion, it is 



so that there is no radiation reaction [cf. eqs. fl25l) or í|26p ]. However, neither the 
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acceleration energy, nor the energy radiated by the charge per unit of proper time, 
^e 2 a l·L a t± , are zero! 

The difference is that in the discrete case this acceleration energy does not exist as 
such, being a term proceeding from the discretized expression for the charged particle 
kinetic energy variation. As we have seen in eq. (l22|) . the Schott term appears when 
the variation of the kinetic energy during the fundamental interval of proper time is 
expanded in powers of r : 

T(r) - r(r - to) ^ _d_ y _ 2 ^ d_ 
To dr 3 c 2 dr 

This is an interesting result, since it was not easy to understand the physical meaning 
of the Schott acceleration energy. With the introduction of the fundamental interval of 
time, as we know, the changes in the kinetic energy are no longer continuous, and the 
Schott term merely expresses, to first order, the variation of the kinetic energy when 
passing from one discrete instant of time to the subsequent one. 

In eqs. 022p and Q25]) . the derivative dT/dr is a point function, forwarding the kinetic 
energy slope at the instant r. And the dissipative term |e 2 a^a M is just a relativistic 
gener alization of the Larmor radiation law: then, if there is acceleration, there is also 
radiation emission. 

For the hyperbolic motion, however, the energy dissipated (because of the accel- 
eration) has just the same magnitude as the energy gain due to the kinetic energy 
increase. We are not forced to resort to 'pre-accelerations' in order to justify the origin 
of such energies. ■·|42j Thus, the present theory provi des a clear picture of the physical 
processes involved in the uniformly accelerated motion of a charged particle. 

3 The Hypothesis of the Chronon in Quantum 
Mechanics 

Let us pass to the main tòpic of the present paper: the chronon in Quantum Me- 
chanics. The speculations about the discreteness of time (on the basis of of possible 
physical evidences) in QM go back to the first dècades of this century, and various 
theories have been proposed developing QM on a space-time lattice (cf., e.g., refs.p£3"]). 
This is not the case with the hypothesis of the chronon, where we do not actually have 
a discretization of the time coordinate. In the twenties, for example, Pokrowski[70j 
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suggested the introduction of a fundamental interval of time, starting from an analysis 
of the shortest wavelengths detected at triat time in còsmic radiation. More recently, 
for instance, Ehrlich[4*4"] proposed a quantization of the elementary particle lifetimes, 
suggesting the value 4.4 x 10~ 24 s for the quantum of time. But a time discretization 
is suggested by the very foundations of QM. There are physical limits that prevent 
the distinction of arbitrarily close successive states in the time evolution of a quantum 
system. Basically, such limitations result from the Heisenberg relations: in such a way 
that, if a discretization is to be introduced in the description of a quantum system, it 
cannot possess a universal value, since those limitations depend on the characteristics 
of the particular system under consideration. In other words, the value of the funda- 
mental interval of time has to change a priori from system to system. All these points 
make the extension of Caldirola's procedure to QM justifiable. 

In the seventies, Caldirola extended the introduction of the chronon to QM, follow- 
ing the same guidelines that had led him to his theory of the electron. So, time is still 
a continuous variable, but the evolution of the system along its world line is discontin- 
uous. As for the electron theory in the non-relativistic limit, one has to substitute the 
corresponding finite-difference expression for the time derivatives; e.g.: 

df (í) f(t)-f(t-At) 

~dT At [ } 

where proper time is now replaced by the local time t. Such a procedure was then 
applied to obtain the finite-difference form of the Schròdinger equation. As for the 
electron case, there are three different ways to perform the discretization, and three 
"Schròdinger equations" can be obtained[4"5]: 



i- [tf (x, t) - V (x, t - r)] = H* (x, t) , (34) 

T 

i^- [* (x, t + r) - * (x, t - t)\ = H* (x, t) , (35) 
Zt 

i- [* (x, t + r) - * (x, t)] = (x, t) , (36) 
r 

which are, respectively, the retarded, symmetric and advanced Schròdinger equations, 
all of them transforming into the (same) continuous equation when the fundamental 
interval of time (that can now be called just r) goes to zero. It can be immediately 
observed that the symmetric equation is of the second order, while the other two are 
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first order equations. As in the continuous case, for a finite-difference equation of order 
n a single and complete solution requires n initial conditions to be specified. 

As the equations are different, the solutions they provi de are also fundamentally 
different. In order to study the properties of such equations there are two bàsic proce- 
dures. For some special cases, they can be solved by one of the various existing methods 
of solving finite-difference equations, or by means of an attempt solution, an ansatz. 
The other way is to find a new Hamiltonian H such that a new continuous Schròdinger 
equation 

ffl ^|ií> = M(x,t) (37) 

reproduces, at the points t = nr, the same results obtained from the discretized equa- 
tions. As has been shown by Casagrande and Montaldi[46j, it is always possible to find 
a continuous generating function which makes it possible to obtain a differential equa- 
tion equivalent to the original finite-difference one, such that in every point of interest 
their solutions are identical. This procedure is useful since it is generally very diflicult 
to work with the finite-difference equations on a qualitative basis. Except for some 
very special cases, they can be only numerically solved. This equivalent Hamiltonian 
H is, however, non-hermitian and it is frequently very diflicult to be obtained. For the 
special case where the Hamiltonian is time independent, the equivalent Hamiltonian is 
quite easy to calculate. For the symmetric equation, e.g., it is given by: 

H =-sm- 1 C-H). (38) 
r n 

As expected, H — > H when r — > 0. Caldirola[47] used the symmetric equation to 
describe the non-radiating electron (bound electron) since for Hamiltonians explicitly 
independent of time its solutions are always of oscillating character: 

# (x, t) = exp (-i- sm^ijÈ)^ f (x) . 

In the classical theory of electron, the symmetric equation also represents a non- 
radiating motion. It provides only an approximate description of the motion without 
taking into account the effects due to the self fields of the electron. However, in the 
quantum theory it plays a fundamental role. In the discrete formalism, it is the only 
way to describe a bound non-radiating particle. 

The solutions of the advanced and retarded equations show completely different 
behaviour. For a hamiltonian explicitly independent of time the solutions have a general 
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form given by 



tt(x,t) 



n 



l-t/r 



/(*) 



and, expanding f(x) in terms of the eigenfunctions of H, 



Hu n (x) = W n u n (x) 



with 



/ ( x ) = Z c « M « ( x ) 



it can be obtained that 



h 



*(x,t)=£c 

n 

In particular, the norm of this solution is given by 



-t/T 



U n (x) . 



|* (x,t)| 2 = Z |c«| 2 exp (-7„í) 



with 



1 

T 



T 



The presence of a damping factor depending critically on the value r of the chronon 
must be remarked. 

This dissipative behaviour originates from the retarded character of the equation. 
The analogy with the electron theory also holds and the retarded equation possesses 
intrinsically dissipative solutions representing a radiating system. The Hamiltonian has 
the same status as in the continuous case: it is an observable since it is a hermitian 
operator and its eigenvectors form a basis of the state space. However, due to the 
damping term, the norm of the state vector is not constant anymore. An opposite 
behaviour is observed for the solutions of the advanced equation, in the sense that they 
increase exponentially. 
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Before going on, let us at least mention that the discretized QM (as well as Caldirola 
et al.'s approach to "QM with friction") can find room within the theories[l8] based on 
the so-calles Lie- admissible algebras. A lot of related work (not covered in the present 
review) can be found e.g. in refs. [i9l \5Ü\ H8] ; see also foT], |52| |53| [51]. 

For a different approach to decaying states, see e.g. [53]. 

3.1 The mass of the muon 

The most impressive achievement due to the introduction of the chronon hypothesis 
in the realm of QM is obtained in the description of a bound electron using the new 
formalism. Bound states are described by the symmetric Schròdinger equation and 
a Hamiltonian that does not depend explicitly on time. A general solution can be 
obtained by using a convenient ansatz: 

* (x, t) = u n (x) exp (-ia n t) , 

n 

where Hu n (x) = E n u n (x) gives the spectrum of eigenvalues of the Hamiltonian. If 
the fundamental interval of time r corresponds to the chronon 8 associated with the 
classical electron, it can be straightforwardly obtained that 

1 . -i / E n 6 \ 

an = T sm {—)■ 

This solution gives rise to an upper limit for the eigenvalues of the Hamiltonian due 
to the condition 



E n 9 



< 1. 



k 

Since #o is fmite, there is a maximum value for the energy of the electron given by 

Elm » 2 -^t „ 105.04 MeV. 
3 e 2 

Now, including the rest energy of the electron, we finally get 

E = £ max + ^ lectron w 105.55 MeV. 

which is very close (an error of 0.1%) to the measured value of the rest mass of the 
muon. Using the equivalent hamiltonian method it is possible to extend the basis of 
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eigenstates out of the critical limit. However, for the eigenvalues above the critical limit 
the corresponding eigenstates are unstable and decay in time: 

* (x, t) = c n u n (x) exp (-ij n t) exp (-k n t) , 

n 

As for the retarded equation, the norm of the state vector is not constant and decays 
exponentially with time for those eigenstates out of the stability range. Caldirola[47j 
interpreted this norm as giving the probability of the existence of the particle in its 
original Hilbert space and associated a mean lifetime with these states. 

The considerations regarding the muon as an excited state of the electron can be 
traced back to the very days of its discovery. Particularly, it has already been observed 
that the ratio between the masses of the two partides is almost exactly 3 /(2a), where 
a is the fine structure constant. [57] It has already been remarked that |a is just the 
coefficient of the radiative reaction term in Dirac's equation for the classical electron. [58J 
Bohm and Weinstein[33] put forward the hypothesis that various kinds of "mesons" 
could be excited states of the electron. Dirac[59] even proposed a specific model for an 
extended electron so as to intèrpret the muon as an excited state of the electron (on 
this point, cf. also refs.[60j). 

Caldirola[56] observed that by means of the Heisenberg uncertainty relations it is 
possible to associate the existence of the muon as an excited state of the electron with 
the introduction of the chronon in the theory of electron. The relation 

At AE > h/2 

imposes limitations in the determination, at a certain instant r, of the energy E as- 
sociated with the infernal motion of the electron. If excited states of the particle 
corresponding to larger vàlues of mass exist, then it is possible only to speak of an 
"electron with rest mass m " when AE < (fi — m )c 2 , where fio is the rest mass of 
the infernal excited state. Such infernal states could only be excited in the presence of 
sufficiently strong interactions. From the uncertainty relation we have that 



2 (fi - m ) c 2 ' 
and, supposing the muon as an excited state, we get 

/ \ 2 r^j ^ hc o 

(fio -m )c = -^rriQC . 

2 e z 
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Thus, it can be finally obtained that 



Ar > -- 



To 

2 : 



3 m c 2 

i.e., that the value of the rest mass of an interacting electron can be taken only inside 
an interval of the proper time larger than half a chronon. So, when we take into account 
two successive states, each one endowed with the same uncertainty Ar, they must then 
be separated by a time interval of at least 2 Ar, which corresponds exactly to the 
chronon r . 



3.2 The mass spectrum of leptons 

In order to obtain the mass of another particle, a possibility to be considered is to take 
the symmetric equation as describing the muon. According to this naïve argument ation, 
the equation also foresees a maximum limit for the energy of the eigenstates of the muon. 
By assuming the equation as successively describing the partides corresponding to these 
màxima, an expression can be set up for the various limit vàlues, given by 



E, 



{<>) 



such that, for 



m c 



3 hc 
2~é 



+ 1 



m c 



3 1 

2a 



+ 1 



(39) 



n = o -> = 0.511 MeV 

n = 1 -> E<Ü = 105.55 MeV 
n = 2 -> = 21801.54 MeV 



(electron) 
(muon) 
(heavy lepton?) 



the masses for the first excited states can be obtained, including a possible heavy lepton 
which, according to the experimental results up to now, does not seem to exist. 

Following a suggestion of Barut [6T] . according to which it should be possible to 
obtain the excited states of the electron from the coupling of its intrinsic magnètic 
moment with its self field, Caldirola |62|, [63] . considering a model of the extended electron 
as a micro-universe, succeeded in evaluating also the mass of the lepton r . 

Caldirola took into account, for the electron, a model of a point-object moving 
around in a 4-dimensional de Sitter micro-universe characterized by 



2<2 2 2 2 

ert —x—y—z 



1J2 



C T, 



■ 
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where tq is the chronon associated with the electron and the radius of the micro-universe 
is given by a = ct . Considering the spectrum of excited states obtained from the naïve 
argumentation above, we find that each excited state determines a characteristic radius 
for the micro-universe. Thus, for each particle, the trajectory of the point-object is 
confined to a spherical shell defined by its characteristic radius and by the characteristic 
radius of its excited state. For the electron, e.g., the point-object moves around, inside 
the spherical shell defined by its corresponding radius and by the one associated with 
its excited state: the muon. Such radii are given by 



2a 



(40) 



According to the model — supposing that the intrinsic energy of the lepton e*™-* is 



given by m 



the lepton moves in its associated micro-universe along a circular 



trajectory with a velocity (5 — to which corresponds an intrinsic magnètic moment 



u (n) 



1 



(41) 

4 7rm (n) c 2· V ) 

Starting from Barut's suggestion, Caldirola obtained for the lepton an extra 
self-energy given by 



E (n, P ) = ( 2 p) 4 m («) c 2 . 

The condition set down on the trajectory of the point-object, so that it remains 
confined to its corresponding spherical shell, is given by 



+ i 

2a 



m Q c 2 , 



and the vàlues attainable by p are: p = for n = 0, and p = 0, 1 for n ^ 0. The 
spectrum of mass is then finally given by 



m (n, P ) = 



1 + (2pY 

Thus, for different vàlues of n and p we have: 



i + (2 P y 



3 1 

2a 



+ 1 



(42) 



n 


P 


m (n) 










0.511 MeV 


electron 


1 





105.55 MeV 


muon 




1 


1794.33 MeV 


tau 
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It must be remarked that the tau appears as an internal excited state of the muon and 
its mass is in fair agreement [64j with the experimental vàlues: m T pd 1784 MeV. The 
difference between these vàlues is less than 1%. This is quite amazing if we consider the 
simplicity of the model. The model foresees the existence of other excited states which 
do not seem to exist. This is to some extent justifiable once the muon is obtained as 
an excited electron and the description of the electron does not anticipate the existence 
of any other state. In order to obtain the lepton tau it was necessary to introduce in 
the formalism the coupling of the intrinsic magnètic moment with the self-field of the 
electron. 

3.3 Feynman path integrals 

The discretized Schròdinger equations can easily be obtained using Feynman's path 
integral approach. This is particularly interesting since it gives a clearer idea of the 
meaning of these equations. According to the hypothesis of a chronon, time is still a 
continuous variable and the introduction of the fundamental interval of time is con- 
nected — let us recali- only with the reaction of the system to the action of a force. It 
is convenient to restrict the derivation to the one-dimensional case, considering a par- 
ticle under the action of a potential V(x,t). While the time coordinate is continuous, 
we assume a discretization of the system (particle) position corresponding to instants 
separated by time intervals r [cf . figure 13. 3j . 

The transition amplitude for a particle going from an initial point (xi,ti) of the 
space-time to a final point (x n ,t n ) is given by the propagator 



In Feynman's approach this transition amplitude is associated with a path integral, 
where the classical action plays a fundamental role. It is convenient to introduce the 
notation 



such that C (x,x) is the classical Lagrangian and S(n,rii) is the classical action. Thus, 
for two consecutive instants of time, the propagator is given by 




(43) 



f, 



n 




(44) 




(45) 
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Figure 1: Discrete steps in the time evolution of the considered system (particle). 



The path integral has the meaning of a sum over all the possible paths traversed by 
the particle and can be written as 



(x n ,t n | Xi,ti) = lim A N J dxjv-i J àxN-2 •• • J dx 2 JJ exp {^S (n,n — (46) 

where A is a normalization factor. In order to obtain the discretized Schròdinger equa- 
tions we have to consider the evolution of a quantum state between two consecutive 
configurations (a; n _i,£ n _i) e (x n ,t n ). The state of the system at t n is: 

+00 

n-l)^ (%n-l, t n -l) dx n _i- (47) 

—00 

On the other side, it follows from the definition of the classical action (eq. Hlj) that 

S (x n , t n ; x n „i, í n _i) = ^ (x n - x n _i) 2 - tV ( Xn ^ n ~ 1 , t n -ij • (48) 
Thus, the state at t n is given by 

+00 

W(x n ,t n )= J exp j— (x„ - x n _i) g »*»*-iJf 

—00 

^(x n _i,í n _i)dx n _i. (49) 

When t 0, for x ra slightly different from x n -i, the integral due to the quadratic 
term is rather small. The contributions are considerable only for x n ~ x n _i. Thus, we 
can make the following approximation 

X n -i = X n + 7] -> dx n „i = df], 

such that 



* (ar„_i, í„_i) = * (z n , í n _i) + ( v 21 " )v+[^y\V 



'd* ( x w ,t w _i) \ + 



2 



By inserting this expression into equation (1491) . supposing that§ 



^The potential is supposed to vary slowly with 
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and taking into account only the terms to the first order in r, we get 



Notwithstanding the fact that exp(— irV(x n , t n )/K) is a function defined only for 
certain well-determined vàlues, it can be expanded in powers of r, around an arbitrary 
position (x n ,t n ). Choosing A = (2ihirr /m) -1 / 2 , such that r — > in the continuous 
limit, we obtain 



2 

(x n , t n -i + t) - * (x n , í n _i) = -j^V (x n , í n _i) * (x n , í n _i) + 



2m cte 2 

By a simple reordering of terms we finally get 



l = + V ( X n, í„-l)| * (X n , í„_i) 

Following this procedure we obtain the advanced finite-difference Schròdinger 
equation which describes a particle performing a one-dimensional motion under the 
potential V(x, t). 

The solutions of the advanced equation show an amplification factor which may 
suggest that the particle absorbs energy from the field described by the Hamiltonian 
in order to evolve in time. In the classical domain the advanced equation is inter- 
preted as describing a positron. However, in the realm of the non-relativistic QM, it 
is more naturally interpreted as representing a system which takes in energy from the 
environment. 

In order to obtain the discrete Schròdinger equation only the terms to the first order 
in t have been taken into account. Since the limit r — > has not been accomplished, 
the equation thus obtained is only an approximation. The meaning of this fact may be 
related to another one we are going to face below in this paper, when considering the 
measurement problem in QM. 

It is rather interesting to remark that in order to obtain the retarded equation we 
have to regard the propagator as acting backward in time. The conventional procedure 
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used in the continuous case always provides us with the advanced equation. Therefore, 
the potential describes a mechanism for transferring energy from a field to the system. 
The retarded equation can be obtained only by assuming an inversion of the time order, 
considering the expression 



which can be rigorously obtained by merely using the closure relation for the eigenstates 
of the position operator and then redefming the propagator in the inverse time order. 
With this expression, it is possible to obtain the retarded Schròdinger equation. The 
symmetric equation can easily be obtained by a similar procedure. 

An interesting characteristic related to these apparently opposed equations is the 
impossibility of obtaining one from the other by a simple time inversion. The time 
order in the propagators must be related to the inclusion, in these propagators, of 
something like the advanced and retarded potentials. Thus, in order to obtain the 
retarded equation we have to take into account effects that act backward in time. 
Considerations like these that led to the derivation of the three discretized equations 
can supply useful guidelines for the comprehension of their meaning. 

3.4 The Schròdinger and Heisenberg pictures 

In discrete QM, as well as in the "continuous" one, the utilization of discretized Heisen- 
berg equations is expected to be preferable for certain types of problems. As it happens 
for the continuous case, the discretized versions of the Schròdinger and Heisenberg 
pictures are also equivalent. However, we show below that the Heisenberg equations 
cannot, in general, be obtained by a direct discretization of the continuous equations. 

First of all, it is convenient to introduce the discrete time evolution operator for the 
symmetric 




(51) 




(52) 



and for the retarded equation, 



Ü(t,t )= 1 + ttH 



ío)/r 



(53) 
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In order to simplify the equations, the following notation will be used throughout 
this Section 



A/(í) 



A fl /(í) 



f(t + T)-f(t-T) 

2t 

, f(t)-f(t-r) 

T 



(54) 



(55) 



For both the operators above it can easily be demonstrated that, if the Hamiltonian 
H is a hermitian operator, the following equations are vàlid: 



AU (Mo) = ^(Mo)#, 



AC/t (í,í ) = --t>t ( Mo )#. 



(56) 



(57) 



In the Heisenberg picture the time evolution is transferred from the state vector to 
the operator representing the observable according to the definition 

i H = & (t,t Q = 0)À s Ü(t,t = 0) (58) 

For the symmetric case, for a given operator A s , the time evolution of the operator 
A ü (t) is given by 



AÍ H (t) = A W (t, t = 0) À S Ü (t, t = 0) 



1 



AA R (t) = — A H ,H 



ih 



(59) 



which has exactly the same form as the equivalent equation for the continuous case. The 
important feature of the time evolution operator which is used to derive the expression 
above is that it is an unitary operator. This is true for the symmetric case. For the 
retarded case, however, this property is not satisfied anymore. Differently from the 
symmetric and continuous cases, the state of the system is also time dependent in the 
retarded Heisenberg picture: 



* H (Í) 



1 + 



r 2 H 2 



h 2 



-(t-* )/T 



(60) 



By using the property À, f (^Aj = 0, it is possible to show that the evolution law 



for the operators in the retarded case is given by: 
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{ 



1 r 



} 



H 



AÀ 



(t) 



À s (t),H s (t) +AÀ S 



(t) 



(61) 



In short, we can conclude that the discrete symmetric case and the continuous case 
are formally very similar and the Heisenberg equation can be obtained by a direct 
discretization of the continuous equation. For the retarded and advanced cases, how- 
ever, this does not hold. In the Appendices we analyse the compatibility between the 
Heisenberg and Schrodinger pictures. 

Let us here mention that a lot of parallel work has been done by Jannussis et 
al.; they, e.g., studied the retarded, dissipative case in the Heisenberg representat ion, 
passing then to study in that picture the (normal or damped) harmònic oscillator: 
see[5T] (cf. also[52l [53]). 

3.5 Time-dependent Hamiltonians 

When we restricted the analysis of the discretized equations to the time independent 
Hamiltonians this was made aiming at simplicity. When the Hamiltonian is explicitly 
time dependent the situation is very similar to the continuous case. It is always difficult 
to work with such Hamiltonians but, as in the continuous case, the theory of small 
perturbations can also be applied. For the symmetric equation, when the Hamiltonian 
is of the form 



such that V is a small perturbation related to H , the resolution method turns out to be 
very similar to the usual one. The solutions are equivalent to the continuous solutions 
followed by an exponentially varying term. It is always possible to solve this kind of 
problem using an appropriate ansatz. 

However, another factor must be considered, related to the existence of a limit 
beyond which H does not have stable eigenstates. For the symmetric equation, the 
equivalent Hamiltonian is given by 



Thus, as previously stressed, beyond the critical value the eigenvalues are not real 
and the operator H is not hermitian anymore. Below that limit, H is a densely defined 
and a self-adjoint operator in the C C L 2 subspace defined by the eigenfunctions of H. 



H = H + V(t), 



(62) 




(63) 
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When the limit value is exceeded the system changes to an excited state and the previous 
state loses physical meaning. In this way, it is convenient to restrict the observables 
to self-adjoint operators that keep invariant the subspace C. The perturbation V is 
assumed to satisfy this requirement. 

In the usual QM it is convenient, in order to deal with time dependent perturbations, 
to work with the interaction representation (Dirac's picture). In this representat ion, 
the evolution of the state is determined by the time dependent potential V(t), while the 
evolution of the observable is determined by the stationary part of the Hamiltonian H . 
In the discrete formalism, the time evolution operator defined for H , in the symmetric 
case, is given by 



Uq (t, t ) = exp 



i(t-t ) 



sm 



(64) 



r \ h 

In the interaction picture the vector state is defined, from the state in the 
Schròdinger picture, as 



*\t)\ = UÍ(t) * s (í ) 



(65) 



where Üç,(t) = Üç,(t, t = 0). On the other hand, the operators are defined as 



A l = Ul(t) A s U (t) 



(66) 



So, it is possible to show that, in the interaction picture, the evolution of the vector 
state is determined by the equation 



ihA^ 1 (x,t) 



ih 

27 



tf 1 (x, t + r) - tf 1 (x, t - r)l = V 1 ^ 1 (x, t) , 



(67) 



which is equivalent to a direct discretization of the continuous equation. For the oper- 
ators we get that 



Ai 1 (í) 



À 1 (t + r) - À 1 (t - r) 



(68) 



2r ih 
which is also equivalent to the continuous equation. 

Thus, for the symmetric case, the discrete interaction picture keeps the same char- 

acteristics of the continuous one for the evolution of the operators and state vectors 

once, obviously, the eigenstates of H remain below the stability limit. We can adopt, 

for the discrete procedure similar to that one commonly used in QM to deal 

with small time dependent perturbations. 
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We consider, in the interaction picture, the same basis of eigenstates associated with 
the stationary Hamiltonian H , given by \n). Then, 

M*)) 1 = £*(*) (n\*(t)) 1 \n) = j:c n (t)\n) 

n n 

is the expansion, over this basis, of the state of the system at a certain instant t. It 
must be observed that the evolution of the state of the system is determined once the 
coefficients c n (t) are known. Using the evolution equation (|67|) it can be obtained that 

ihA (n | * (í)) 1 = J2 H V* \ m ) ( m I * (*))*■ 

m 

Using the evolution operator to rewrite the perturbation V in the Schròdinger picture 
we get that 



ÍhAc n (t) = C m (í) V nm (t) exp (ÍUJ nm t) 



(69) 



such that 



T 



. _ x (rE n \ . _ x {tE„ 
sm — — sm 



h )~~" V h 

and we obtain the evolution equation for the coefficients c n (t), the solution of which 
gives the time evolution of the system. 

As in the usual QM, it is also possible to work with the interaction picture evolution 
operator, U l (t,t ), which is defined as 

such that (1671) can be written as 



ihAÜ\t,t ) =V\t)Ü\t,t ). (70) 

The operator Ü l (t,t ) has to satisfy the initial condition U l {t, to) = 0. Given this 
condition, we have for the finite-difference equation above the solution 



Ü\t,t 



o) = exp 



-z(í-ío) • -i frV\ty 
sm 1 



r V h 

Differently from the continuous case, where the approximate evolution operator 
turns out to be an infmite Dyson series, a well determined expression is obtained. The 
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solution of the problem is obtained by correlating the elements of the matrix associated 
with such operator to the evolution coeficients c n (t). 

In general, the finite-difference equations are harder to be analytically solved than 
the equivalent differential equations. In particular, such difficulty is far more stressed 
for the system of equations obtained from the formalism above. 

An alternative approach is to use the equivalent Hamiltonians. [56] Once the equiv- 
alent Hamiltonian is found the procedure is exactly the same as for the continuous 
theory. If the perturbation term V is small the equivalent Hamiltonian can be written 
as 



H 



V{t) = È + V(t). 



In the interaction picture, the state of the system is now defined as 



exp %- 



and the operators are given by 



ít / .HqÍ \ jc 

A = exp % A exp 

\ h J 



The states ( J71j) evolve according to the equation 





lh di 



m\t)) = v 1 v\t) 



(71) 



(72) 



(73) 



where V 1 is obtained according to definition (1721) . 

Now, small time dependent perturbations can be dealt with by taking into account 
the time evolution operator defined by 



y\t)) = Ü\t,t ) ¥(t )). 
According to the evolution law (1731) we have 



(74) 



zh-Ü\t,t )=V\t)Ü\t,t ) 



(75) 



Thus, once given that ^(ío^o) = 1 it turns out that the time evolution operator is 
given by 



Ü I {t,t Q ) = l-^ f V\t')Ü\t\t Q )<\t' 
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or 



Í/ I (t,í ) = l + E (-*) / d W dí 2 ···/ àt n V 1 {t 1 )V\t 2 )···V\t n ), 

~y\ TiJ Jta Jto Jt 

where the evolution operator is obtained in terms of a Dyson series. 

In order to draw a parallel between the elements of the matrix of the evolution 
operator and the evolution coefficients c n (t) obtained from the continuous equation 
equivalent to (|69|) . we have to use the basis of eigenstates of the stationary Hamiltonian 
Ho- If the initial state of the system is an eigenstate \m) of that operator then, at a 
subsequent time, we have 

c n (t) = (n\Ü\t,t )\m) . 

The method of the equivalent Hamiltonian is simpler to use, since it takes full 
advantage of the continuous formalism. 

4 Some Applications of the Discretized 
Quantum Equations 

Turning back to more general qüestions, it is interesting to analyse the physical 
consequences resulting from the introduction of the fundamental interval of time in 
QM. In this Section we apply the discretized equations to some typical problems. 

4.1 The simple harmònic oscillator 

The Hamiltonian that describes a simple harmònic oscillator does not depend explicitly 
on time. The introduction of the discretization in the time coordinate does not affect the 
outputs obtained from the continuous equation for the spatial branch of the solution. 
This is always the case when the potential does not have an explicit time dependence. 
For potentials like this, the solutions of the discrete equations are always formally 
identical, with changes in the numerical vàlues which depend on the eigenvalues of the 
Hamiltonian considered and on the value of the chronon associated with the system 
described. We have the same spectrum of eigenvalues and the same basis of eigenstates 
but with the time evolution given by a different expression. 

For the simple harmònic oscillator, the Hamiltonian is given by 
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H 



1 mu 2 - 2 
P 2 H — X 2 



2m 2 
to which the eigenvalue equation corresponds: 



(76) 



H \u r . 



E n \u n ) 



(77) 



so that E n gives the energy eigenvalue spectrum of the oscillator. 

As mentioned previously, since this Hamiltonian does not depends explicitly on time, 
there is always an upper limit for the possible vàlues of its energy eigenvalues. In the 
basis of eigenfunctions of H a general state of the oscillator can be written as 



W)> = E c «(°) K) ex P 



.t . 

—i — sm 

T 



E n r 

h 



with c n (0) = (u n \^f(t = 0)). Naturally, when r — > 0, the solution above recovers the 
continuous expression with its time dependency given by exp ( ~^ nt ) • Therefore, there 
is only a small phase difference between the two expressions. For the mean value of an 
arbitrary observable, 



*(f)l4*(*)) = EE4(oK(ok 

m=0 n=0 



exp 



exp 



( E m E n) m 3 



+ 0(r 4 ) , 



with A mri = (u m \A\u n J, we obtain an additional phase term which implies a small 
deviation of the resulting freqüències when compared to the Bohr freqüències of the 
harmònic oscillator. To first approximation, this deviation is given by the term de- 
pending on r 2 in the expression above. 

It must be emphasized that the restrictions imposed on the spectrum of eigenvalues 
of H break the basis of eigenvectors: the number of eigenvectors becomes finite and 
does not constitute a complete set and, therefore, does not form a basis anymore. For 
eigenstates beyond the upper limit the states are unstable and decay exponentially with 
time. 

For a time independent Hamiltonian, the retarded equation always furnishes damped 
solutions characteristic of radiating systems. In this case there is neither stationary 
solutions nor upper limit for the energy eigenvalues. The larger the eigenvalue the 
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larger the damping factor and more quickly its contribution to the state of the system 
tends to zero. If we write the state of the oscillator as 



!*(*)> = E c «(°) K) 

n 

which has a norm decaying according to 



1 + TTE n 



<*(í)|*(í)) = E c "(0)K> 



T 2 E 2 
1 + — 

+ h 2 



we have for an arbitrary observable that [with (A(t)) = (A) (i)]: 



(78) 



(M*)) = EE c m(0)c„(0)A mn exp 

m n 

or, to the first order in r, 



— In 



T %T 

1 + -TïE n E m — — (E m — E n ) 
n n 



W)) = EE^(0)c n (0)A mn exp 



(E m — E n ) t 



cxp 



i El 



El 



T 

2h 2 



so that, beside the Bohr freqüències defming the emission and absorption freqüències 
of the oscillator, we obtain a damping term which causes the average value of the 
observable — which is explicitly independent of time — to tend to zero with time. A 
cursory analysis shows that even for very small eigenvalues, smaller then 1.0 eV, the 
damping factor is large, so that the decay of the average vàlues is very fast. The 
damping factor of the norm in equation (178"]) can be evaluated, and its behaviour can 
be seen in figure 14.11 



4.2 Pree particle 

For a free particle, an electron for example, the general solution of the symmetric 
equation (135]) can be obtained, in the coordinate representat ion, using as an ansatz the 
solution for the continuous case. Thus, a spectrum of eigenfunctions (plane waves) is 
obtained given by 

%(x,t) = (27thy 3/2 exp (-ia(| p |)í + zÍP-^y 
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t(s) 
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Figure 2: Typical behaviour of the damping factor associated with different energy 
eigenvalues (retarded case). 



Inserting this expression into the symmetric equation, we obtain for the frequency 
a(|p|) that 



When r — > 0, a(|p|)7i coincides with the energy of the particle. As has been observed 
for the bound particle, here we also have an upper limit for the spectrum of eigenvalues. 
Thus the upper limit for the possible vàlues of momentum is given by 



V < PMax = y — f~ = 10 MeV/c (80) 

for the electron. In other words, there is a limit beyond which the freqüències cease to 
be real. 

As in the continuous case, the state of the particle is described by a superposition 
of the eigenstates and can be written as 

1 /' , , _ .. ( ,_/, , ,(P" X A 



^( x ' *) = (t? n] 3/2 J d V(p) exp ( -ia(|p|)í + i 



The coefficients c(p) are determined from the initial condition \l/(x, 0) = \I/o(x). 
From the expression for a, it can be observed that beyond a certain value of p the 



expression loses meaning. When p > y^ 1 , the complete solution will be defined only 
if c(p). From the stationary phase condition we have that 



"" , · /l " (t) 



2 / ^4 



and, supposing that c(p) corresponds to a distribution of probabilities with a peak at 
p = Po , then the wave packet will move in the direction p with uniform velocity 



m 



1 - 



T\ I p 



-1/2 



which coincides with the group velocity of the packet. It can be promptly observed that 
when p reaches its maximum value permitted, the velocity diverges: v — * oo. Thus, 
the introduction of a fundamental interval of time does not bring in any restriction to 
the velocity of the particle, although it results in a limit for the canonical momentum 
of the eigenfunctions. Starting from the condition of stationary phase it is possible to 
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redefine the momentum associated with the particle, so that this new momentum does 
not suffer any restriction at all. So, one can conclude that the existence of free electrons 
with any energy is possible, differently from what happens for the bound electron. 

For p > pMax the frequency «(|p|) fails to be real and its dependence on p is shown 
in figure H~2l An analysis of equation (1751) shows that if a(|p|) is complex then, for 
P < PMax, the imaginary component is null and the real part is given by expression 
(EH])- Whenp>p M ax, then 



Re (oí(jp)) 



7Ï 

2? 



Im (a{p)) = ln 



T 



rp z 



2m h 



+ 



rp z 



\ \2m h / 



1 



with the real part being a constant and the imaginary one tending logarithmically to 
— oo. Using the expressions above we can observe that, for p > pMax, the eigenstates 
become unstable, with a time dependent decay term. When we look for an equiva- 
lent Hamiltonian H that, for the continuous Schrodinger equation, supplies equivalent 
outputs, we have that this is possible only if H is a non-hermitian operator. It is 
straightforward to see that this is the case for H = H\ + ÍH2, with H\ and H2 hermi- 
tian and such that Hi \p) = H Re (a(p)) \p) and H 2 \p) — h Im \p) . 

For the retarded equation, using the same ansatz of the symmetric case, the damping 
factor appears for every value of p. There is no limitation on the vàlues of p but, when 
p — > 00, the real part of a(|p|) tends to the same limit value observed for the symmetric 
case. Figure I4T21 illustrates the behaviour of the components of a(|p|). The general 
expression for an eigenfunction is found to be 



\l/p(x, t) oc exp 



ipx 



it 1 p T 

— tan - 

r 2mn 



exp 



t 

2r 



ln 



1 



ÍÉL· 

\ 2mh 



Performing a Taylor expansion and keeping only the terms to the first order in r we 
obtain the continuous solution multiplied by a damping factor: 



UJ 2 Tt 



1) 



T / \ í W X ■ \ 

W p (x, t) oc exp I — xut j exp I - 

where ui = p 2 /2mh is the frequency obtained for the continuous case. 

It must be remarked that the damping term depends only on the Hamiltonian, 
through the frequency u, and on the chronon associated with the particle. As the 
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Figure 3: Real and imaginary components of a(|p|) obtained for the symmetric equation 
compared to the continuous case. 
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Figure 4: Real and imaginary components of a(|p|) obtained for the retarded equation 
compared with the continuous case. 



latter is constant for a given particle, that term shows that for very high freqüències the 
solutions decay quite fast and, as the system evolves, the decay for smaller freqüències 
also comes true. 

In figure H~2l it can be observed that the inflection point, delimiting the region of the 
spectrum where the decay is faster, moves in the direction of smaller freqüències as time 
goes by. The consequence of this decay is the narrowing of the frequency bandwidth 
which is relevant for the wave packet describing the particle. This is an echo of the 
continuous decrease of the energy. As in the symmetric case, obtaining an equivalent 
Hamiltonian is possible only if non-hermitian operators are considered. 

At this point, it is worthwhile to reconsider the question of the physical meaning 
of the three discretized Schròdinger equations. Apparently, the choice of the equation 
to be used in a particular situation is determined by the boundary conditions, by the 
restrictions they impose on the system. The symmetric equation is used for special 
situations for which the system neither emits nor absorbs radiation, or does it in a 
perfectly "balanced" way: This is the case for the electrons in their 'atòmic orbits'. So, 
the particle is stable until a certain energy limit, beyond which the behaviour of the 
states is similar to that of the retarded solutions. For energies far below that limit, the 
particle behaves almost identically to the continuous case, only that the new freqüències 
associated with each wave function differ from the continuous freqüències by a factor 
of order r 2 . The probability that a particle is found with energy larger than the limit 
value decreases exponentially with time. For the bound electron, the limit is that one 
equivalent to the rest mass of the muon. If a parallel with the classical approach is vàlid, 
the symmetric equation describes an isolated system, which does not exchange energy 
with the surrounding environment; or a situation of perfect thermodynamic equilibrium, 
in which a perfect balance between absorbed and dissipated energies is verified. For the 
classical theory of the electron the symmetric equation is only an approximation which 
ignores the radiation reaction effects. In QM, however, the existence of non-radiating 
states are related to the very essence of the theory. The symmetric equation shows 
that, below the critical limit, the states are physically identical to the outputs from the 
continuous theory: they are non-radiating states. 

The retarded equation represents a system which somehow loses energy for the 
environment. The mechanism of such energy dissipation is related to the hamiltonian 
of the system but also to some property of the environment — even the vacuum — as it 
can be inferred from the description of the free particle. From the solutions obtained 
it is now observed that time has a well-defined direction of flux and that the frequency 
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Figure 5: Displacement with time of the inflection point. 



composition of the wave packet associated with the particle depends on the instant of 
time considered. It is clear that it is always possible to normalize the state in a certain 
instant and consider it as being an initial state. This is permitted by the formalism. 
However, in a strictly rigorous description, the frequency spectrum corresponds to a 
specific instant of time which took place after the emission. This is an aspect that can 
be interesting from the point of view of possible experimental verifications. 

4.3 The discretized Klein— Gordon equation 

Another interesting application is the description of a free scalar particle — a scalar 
or zero spin photon, — using a finite-difference form of the Klein-Gordon equation for 
massless partides. 

In the symmetric form the equation is written as 

= — > -i ^ - V 2 ^(í) = 0. (82) 

Using a convenient ansatz we obtain, for this equation, in the coordinate representat ion, 
that 

^ k {x,t) = Aexp (-277- cos _1 (l - 2c 2 r 2 k 2 )j exp (ikx), 
which can be written as 

ty p (x,t) = Aexp ^-i^-cos _1 (l - 2c 2 T 2 E 2 /h 2 )^j exp (ipx/h), 

since E = p 2 c 2 and p = Tik. Expanding the time exponential in powers of r, we find that, 
to the second order in r, a solution which is very similar to the continuous expression: 



%(x,t) = Aexp (E't-px)^, 



with 

E' w E 1 



6h 2 ) ■ 

A difference of the order of r 2 is observed between the energy vàlues of the photons 
in the continuous and discrete approaches. The general solution is given by a linear 
combination of the eigenfunctions found. A priori, the value of the chronon for the 
particle is not known. The time dependent exponential term in the expressions above 
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leads to an upper limit for the allowed energy, which is given by E < h/r. We could 
suppose that the value of the chronon for this 'photon' is of about the fundamental 
time interval of the electromagnètic interactions, around 10~ 9 s , resulting in a critical 
value of approximately 6.6 keV, which is a very low limit. A smaller chronon should 
increase this limit but, if there is any generality in the classical expression obtained for 
the electron, we should expect a larger value for this massless particle. 

Whether, instead of a 'photon' we consider a scalar 'neutrino' , taking for the value 
of the chronon r ~ 1CT 13 s — typical time for the weak decay — , the limit for the energy 
associated with the eigenfunctions is now approximately 0.007 eV. This means that in 
the composition of the wave packet describing this particle the only contribution comes 
from eigenfunctions, the energy of which is below that limit. 

The eigenfunctions obtained for the Hamiltonian considered are "plane waves" So- 
lutions. The dependence of these solutions on energy and time is displayed in figures 
14.31 and 14.31 For smaller vàlues of r the decay of the modes with energy above the 
maximum is faster. 

Apparently, it seems to be possible to determine a limiting value for the chronon 
starting from the uncertainty relations. This could be obtained, when describing parti- 
des, using the expression 

h 

2m c 2 

that provides for the electron a maximum limit given by 6.4 x 10~ 22 s. However, this 
value is two degrees of magnitude larger than the classical value of the chronon for the 
electron, which is a considerable difference. For a complex system it is also possible to 
make use of this relation, as we are going to mention ahead. 

We also have to consider the conditions a 'photon' must be supplied with, in order 
to be described by the symmetric equation. [For the electron, it seems clear that not 
irradiating in a bound state — which is imposed by QM — implies the adoption of the 
symmetrical equation]. For the photon, as for a free particle, when using the retarded 
form of the Klein-Gordon equation, a solution is also obtained wherein the highest 
freqüències decay faster than the lowest ones. There is always a tendency in the sense 
that the lowest freqüències prevail. If we are allowed to assign a physical meaning 
to such a discretized Klein-Gordon equation, we are also allowed to think that, the 
farther the light source, the more the spectrum of the emitted light will be shifted for 
the largest wavelengths, even if the source is at rest with respect to to the observer. 
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Figure 6: Solution of the discretized Klein-Gordon equation, when the energy is smaller 
than the critical limit, depicted for different vàlues of energy and time. 




Figure 7: Solution of the discretized Klein-Gordon equation when the energy is larger 
than the critical limit, depicted for different vàlues of energy and time. In this case the 



Thus, we could obtain a red shíft effect as a consequence of the introduction of the 
chronon which could be used in the construction of a 'tired-light' theory. 

Finally, we have to keep in mind that the discretization considered for the Klein- 
Gordon equation does not follow exactly the same procedure which led to the discretized 
Schròdinger equation, since it is a relativistic invariant equation. We did not change 
the proper time but the time coordinate itself into the discretized form. We considered 
a discretized version of the hamiltonian operator by applying the transformations: 



P 




H — ► ihA, 
with A as defined in Subsection 13.41 on the hamiltonian of a relativistic free particle, 

H = \Jp 2 c 2 + m 2 c 4 

as usual in the continuous case. 



4.4 Time evolution of the position and momentum 
Operators: The harmònic oscillator 

It is possible to apply the discretized equations to determine the time evolution of the 
position and momentum operators, which is rather interesting for the description of 
the simple harmònic oscillator. In order to do that, we use the discretized form of 
the Heisenberg equations which, in the symmetric case, can be obtained by a direct 
discretization of the continuous equation. Starting from this equation we determine the 
coupled Heisenberg equations for the two operators: 

p(í + r) - p(í - r) 2 

= —mu x(í), (83) 



2r 

x(í + r) - x(í - r) 1 



p(í). (84) 



2r m 

Such coupled equations yield two finite-difference equations of second order, the 
general solutions of which are easily obtained. The most immediate way to determine 
the evolution of these operators is to use the creation and annihilation operators. Keep- 
ing the Heisenberg equation and remembering that for the harmònic oscillator we have 
H = üü (À^À + |J, we get for the symmetric case: 
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À{t + r) - À{t - r) 



2r 



-iuM(t), 



55) 



such that 



À^t + r) -At(í -r) 



2r 



it(í), 



56) 



A(í) = i(0) exp ( -z- sin" 1 (wr 



(87) 



i f (í) = i f (0) exp fi- sin" 1 (cur 

where we used the fact that, for t = 0, the Heisenberg and Schròdinger pictures are 
equivalent: À(t = 0) = A = i(0) and À^(t = 0) = i 1 " = At(0), with i and i f 
independent of time. In order to obtain these equations we considered that, for the 
non-relativistic case, there is neither creation nor annihilation of partides, such that we 
can impose restrictions on the freqüències in the phase term of the operators. For 
the creation operators, e.g., the terms with negative freqüències — associated with 
antiparticles — are discarded. 

We can observe that the Number, as well as the hamiltonian operator are not 
altered: 



N = À\t)À{t) = i f (0)i(0), 



H = huj[N + 



i) = nu (it(o)i(o) + \ 

Thus, starting from these operators, we obtain for the symmetric case: 



x(í) = x(0) cos 



t . 
— sm 

r 



, P(0) • 
H sm 



— sin 1 (lot) 
t 



p(t) = p(0) cos 



■ sin 1 (cüt) 



t 



mwx(0) sin 



sin 1 (lut) 



ÍT 



which differ from the continuous case since the frequency u here is replaced by a new 
frequency - sin -1 {lot) that, for r — > 0, tends to the continuous one. Also, there is now 
an upper limit for the possible oscillation freqüències given by uo < 1/r. Above this 
frequency the motion becomes unstable, as it can be observed in figure 14.41 
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Figure 8: Phase space of the harmònic oscillator when uj > - and in the discrete case, 
with time intervals múltiples of r: in case (a), time is regarded as intrinsically discrete, 
so that in the picture only the points wher61he lines touch one another are meaningful, 
while in case (b) time is regarded as intrinsically continuous. In the actually continuous 
case no modification is expected with respect to the ordinary case, (c), under the present 
hypotheses. 



The existence of a maximum limit for the frequency is equivalent to an upper limit 
for the energy eigenvalues given by E n = (n+ \)Tno < h/r, which is equal to the upper 
limit obtained using Schròdinger' s picture. Since 

- sin^ 1 (üüt) ^uj + ^uj 3 t 2 + 0(r 4 ), 
r 3! 

the difference expected in the behaviour of the oscillator with respect to the continuous 
solution is quite small. If we take, for example, the vibration frequency of the hydrogen 
molecule {H 2 ), we have that ui ~ 10 14 Hz, while the term of the second order in r is 
smaller than 10~ 3 Hz (if the analogy with the classical theory is vàlid, the chronon is 
expected to be smaller for more massive systems). In terms of average vàlues we have 
that, for the position operator 

(x(í)) = (x(í)) cont + ^t(p(í)) , 

àlm 

in which the term of order r 2 is expected to be considerably smaller than the mean 
value for the continuous case. At this point, naturally, the mean vàlues are determined 
taking for the system a state made up of a superposition of stationary states. For the 
stationary states \u n ) themselves the mean vàlues of x and p are zero. 

For the retarded case the solutions can be obtained using the time evolution oper- 
ators for the Heisenberg equation (Appendix |A}. As expected, decaying terms come 
out. The creation and annihilation operators obtained for this case are then given by 



A(t) = A(0) 1 + íujt + t 2 uj 2 Ç T ta A(Q) exp {—íujt) exp 



i f (í) = i f (0) [l - íujt + rV^l T w i f (0) exp {íujt) exp 



UJ 2 Tt 



with ^ being a real positive factor. The relation (À^Y = À continues to be vàlid but the 
Number operator and, consequently, the Hamiltonian, is now not constant anymore: 



H(t) = À\0)À(0) \(l+u 2 T 2 tf + 



2 2 
ÜJ T 



H(0) exp 



-2 £ + - uj'rt 



H(t) = hu { i f (0)i(0) 



,2^2 
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Taking into account the terms to the second order in r, we get that the oscillation 
freqüències also decays with time. These results are consistent with the fact that the 
system is emitting radiation, with the conseqüent reduction of its total energy. How- 
ever, it is remarkable that the energy of the quanta associated with the creation and 
annihilation operators is not constant, even with a very tiny variation rate. In the 
same way, when we calculate the position and momentum operators a damping factor 
is obtained. Figure EOl· a shows the strange damping factor associated with the Number 
operator. It can be observed that this damping occurs within a period of time which 
is characteristic for each frequency, being slower and postponed for lower freqüències. 
Figure FOl -b shows the dampening of the oscillations as described by the retarded equa- 
tion. Once the expressions for the position and momentum operators are determined, 
we obtain that, to first order in r, 



(x(í)) = i x(0) cos (ut) + sin (u>r) 1 exp 

meu 



(x(í)) = (x(í)) cont exp 



Taking into account the higher order terms, we can observe a small variation in the 
oscillation frequency just as observed in the symmetrical case. The introduction of time 
independent perturbations does not cause any additional variations apart from those 
found even in the continuous case. It must be pointed out that the results obtained 
with the procedure above are in agreement with those obtained following Schròdinger's 
picture. 

4.5 Hydrogen atom 

The hydrogen atom is basically a system made up of two partides attracting each other 
through coulombian force, which is therefore inversely proportional to the square of the 
distance between them. The bàsic Hamiltonian is given by 

P 2 e 2 

h ° = 2-»-r' (89 > 

and is composed of the kinetic energy of the atom in the centre-of-mass frame, and of 
the coulombian electrostatic potential (/i is the reduced mass of the system electron- 
proton). A more complete description is obtained by adding correction terms (fine 
structure) to the Hamiltonian, including relativistic effects such as the variation of the 
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Figure 9: (a) Damping factors associated with the Number, operator calculated for a 
few freqüències, (b) Damping of the oscillations for the harmònic oscillator described 



electron mass with velocity and the coupling of the intrinsic magnètic moment of the 
electron with the magnètic field due to its orbit (spin-orbit coupling). Besides, there are 
also the hyperfine (hf) corrections which appear due to the interaction of the electron 
with the intrinsic magnètic moment of the proton and, finally, the Lamb shíft, due to 
the interaction of the electron with the fluctuations of the quantized electromagnètic 
field. The Hamiltonian can finally be written as [65] 

P 4 . . 

H Y = m e c 2 + H — — L • S + H M + # Lamb . (90) 

micrHr 

The introduction of the magnètic moment of the nucleus through the hyperfine 
correction causes the total angular momentum to be F = J + 1. The Hamiltonian does 
not depend explicitly on time such that, for the symmetric Schròdinger equation 

i— [* ( x , t + r) - * (x, t - r)] = (x, t) , (91) 
Zt 

we obtain, using the separation of variables, the following uncoupled equations: 

#i$(x) = £$(x) 

i^[T(t + T)-T(t-T)] = H 1 T(t), 



with the general solution 



* (x, t) = $ (x) exp 



- l r sm hr) 



(92) 



The difference related to the continuous case appears only in those aspects involv- 
ing the time evolution of the states. Since the Hamiltonian is time independent, its 
eigenvalues are exactly the same as those obtained in the continuous case [65] 

2 1 2 2 m e c 2 ( n 3\ 4 
£(n,i) ~ m c - ^m e c a - \J^T_ ~ l) a + E ™ + 

A situation where a difference between the two cases can appear is when taking into 
account the probabilities of transition between the eigenstates for an atom submitted to 
a time dependent potential. In the discrete approach it is possible to use the method of 
the equivalent Hamiltonian in order to obtain the transition probabilities. As mentioned 
previously (Subsection l3.5p . the problem is treated using the conventional approximate 
methods for time dependent perturbations. 
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If we consider, for example, the non-relativistic interaction of an atom with an 
electromagnètic field described by the vector potential A(x, í), we have for the low 
intensity limit, in the Coulomb gauge, the Hamiltonian 



H(t) = H l - V(t) = Hi- — À (R, t] ■ P, 

m c c v 



(93) 



where the potential term is taken as being the perturbation. If we consider that the 
potential describes a monochromatic field of a plane wave, then 



A (x, t) = A è 



/ . n • x . \ 

exp lüü mot + exp 



n • x 

-iíü h iiüt 



(94) 



c / V c 
where è is the linear polarization of the field and n is the propagation direction. The 
term depending on (— iujt) corresponds to the absorption of a quantum of radiation Tiuj 
and the (iut) term to stimulated emission. Let us assume that the system is initially in 
an eigenstate |<E>j) of the time independent Hamiltonian. Keeping only the perturbations 
to the first order in V(t), we obtain that 

= ~ f exp(iüü ni t')V ni (t')dt' , 
h Jo 

where u ni in the discrete case is given by 



1 

T 



sm 



— sm 



h ) V h 

Working with the absorption term, we get by contrast that 



ieA 
mlch 



$ n |e™/ c (e • p)|^) J* exp [i(u ni - üü)t']dt' 



Thus, the probability of transition from the initial state |<E>j) to the final state |$/) is 
given by 



Pfi(t) = cf\t) 



2 e 2 \A \ 



$/ | e -n.x/c (é . p) | $ M 2 Ç exp _ w)í ,] df / 

J 



or 



Pfi(t) 



4e 2 \A Q \ 2 
m 2 c 2 h 2 



$ f \e lujh - x/c (è-p)\<S> t 



2 sin 2 [(uüfi - u)t/2} 
(u fi - uj) 2 
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so that the determination of the matrix elements of the spacial term, using the elèctric 
dipole approximation, provides the selection rules for the transitions. What is remark- 
able in this expression is the presence of a resonance showing a larger probability for 
the transition when 



1 

LO = Lüfi 



T 



sm — — — sm 1 



This expression is formally different from the one obtained for the continuous ap- 
proach. When we expand this expression in powers of r, we get that 



- - + ^f^- (96) 



E t -Ei lE}- Ef 

li + 6 h 3 

The first term supplies the Bohr freqüències as in the continuous case; the second, the 
deviation in the freqüències caused by the introduction of the time discretization: 

1E 3 -E 3 2 

Atüfi = J „ r 2 . 

} 6 k 3 

If we consider the chronon of the classical electron, r rs 6.26 x 10~ 24 s, it is possible 
to estimate the deviation in the frequency due to the time discretization. Then, for the 
hydrogen atom, 



Acjfi « 2.289 x KT 2 (e 3 - E?) . 



If we take into account, for example, the transitions corresponding to the first lines of 
the series of Lyman and Balmer, i.e., of the non-disturbed states n = n- x — * n = rif, we 
have 



ni 


n f 


AE (eV) 


v (Hz) 


Az/ D (Hz) 


1 


2 


10.2 


2.465 x 10 15 


~ 10 


1 


3 


12.1 


2.922 x 10 14 


~ 10 


1 


4 


12.75 


3.082 x 10 14 


~ 10 


2 


3 


1.89 


4.566 x 10 14 


< 1 



where AE is the difference of energy between the states, v is the frequency of the photon 
emitted in the transition and Az/d is the frequency deviation due to the discretization. 
Such deviation is always very tiny. We must remember that the hyperfine corrections 
and those due to the Lamb shift are of order of a Gigahertz. For the transition n = 
1 — > n — 2, e.g., the correction due to the Lamb shift is approximately 1.06 GHz. 
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Larger deviations caused by the discretization occurs for mono-electronic àtoms with 
larger atòmic numbers. For the first transition the deviation is of approximately 90 Hz 
for the 2 He, 1.1 kHz for the 3 Li, and 420 kHz for the 6 C. However, these deviations 
are still quite smaller than that one due to the Lamb shift. That is also the case for 
the muonic àtoms. For a muonic atom with a proton as nucleus, using for the chronon 
a value derived from the classical expression for the electron (r M = 3.03 x 10~ 26 s) the 
deviation is of about 1.4 kHz for the transition n = 1 — > n = 2. For that transition the 
frequency of the emitted radiation is approximately 4.58 x 10 17 Hz. 

For the retarded equation, a difference with respect to the symmetrical case is present 
in the time evolution of the states. The procedure is identical to the one used above 
and the general solution is now given by 



*(x,í) = $(x) 

so that the transitions now occur with freqüències given by 



TEr" T 
1 + i T 



—ih 

UJ = UJfi 



T 



irEf\ , / ítEí 
ln 1 + - 1 -In 1 + 



n ! -y n jy (97) 

As results from the characteristics of the retarded equation, this is a complex fre 
quency. The real component of such frequency can be approximated by 



Ef-Ej lEf-E?_, 2 

n 3 h 3 



T 



where the first term is the expression for the continuous case. For the particular tran- 
sition n = 1 — > n = 2 we have that the deviation due to the discretization is of about 
18 Hz. 

The imaginary component, on the other hand, can be approximated by 

_ . > i Ef - E? 
Im(wfi) w - 2 r. 

In the expression for the probability of transition we have the module of an integral 
involving the time dependency of the general solution. In this case, the characteristic 
damping causes the probability to tend to a fixed, non zero value. An example of such 
behaviour is shown in figure 14.51 which shows the variation of the time dependent term 
between an initial instant to — and some hundred chronons later. In order to observe 
the decay of the amplitude factor we have used a larger value for the chronon, of about 
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10~ 18 s. When the chronon is of the order of the one we have been considering for the 
electron the decay is slower. 

As it can be observed, the effect of the time discretization on the emission spectrum 
of the hydrogen is extremely small. Using the expressions obtained above we can 
estimate that, in order that the effect of the time discretization is of the same order of 
the Lamb shift, the chronon associated with the electron should be of about 10~ 18 s, 
far above the classical value (but close to the typical interval of the electromagnètic 
interactions). In any case, it should be remembered that the Lamb shift measurements 
do not seem to be in full agreement |70j with quantum electrodynamics. 



Concluding this Subsection, it is worthwhile to remark that, for a time independent 
Hamiltonian, the outputs obtained in the discrete formalism using the symmetric equa- 
tion are very similar to those from the continuous case. For such Hamiltonians, the 
effect of the discretization appears basically in the freqüències associated with the time 
dependent term of the wave function. As already observed, the difference in the time 
dependency is of the kind 



exp [-iE n (t - t )/h] 



exp 



-.sin" 1 (^)(t-t )/r 



The discretization causes a change in the phase of the eigenstate, which can be 
quite large. The eigenfunctions individually describe stationary states, so that the 
time evolution appears when we have a linear combination of such functions, this way 
describing the state of the system. This state evolves according to 



= E c «(°) ex P 



-i sm 



■'^)(*-*.)/r 



V n 



considering that H\(j> n ) = E n \(p n ) is the eigenvalue equation associated with the Hamil- 
tonian. 

When the stationary states of a particle under, e.g., one-dimensional squared poten- 
tials are studied, the same reflection and transmission coefficients and the same tunnel 
effect are obtained, since they are calculated starting from the stationary states. When 
we consider a linear superposition of these stationary states, building a wave packet, 
the time dependent terms have to be taken into account, resulting in some differences 
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Figure 10: Behaviour of the time dependent component of the transition probability. 



with respect to the continuous case. Some attempts have been carried out in order to 
find out significant measurable differences between the two formalisms[.66] |67] but no 
encouraging case has been found yet. 



5 Density Operators and the "Coarse Graining" 
Hypothesis 

5.1 The "coarse graining" hypothesis 

First of all, it is convenient to present a brief review of some tòpics related to the 
introduction of the coarse grained description of a physical system. This hypothesis 
is then going to be used to obtain a discretized form of the Liouville equation, which 
represents the evolution law of the density operators in the usual QM. 

An important point to be remarked is that the introduction of a fundamental inter- 
val of time is perfectly compatible with a coarse grained description. The bàsic premise 
of such description, in statistical physics, is the impossibility of a precise determination 
of the position and momentum of each particle forming the system, in a certain instant 
of time. Let us consider, for the sake of simplicity, a system composed by N similar 
pointlike partides, each of them with three degrees of freedom described by the coordi- 
nates {qi, q 2 , q^) ■ We can associate with this ensemble of partides an individual phase 
space (named /x-space) defined by the six coordinates (çi, q 2 , q3',Pi,P2,P3) so that the 
system as a whole is represented by a crowd of points in this space. 

Since the macroscopic observation is unable to precisely determine the six coordi- 
nates for each particle, let us assume that it is possible only to know if a given particle 
has its coordinates inside the intervals (qi + dqi) and (pi + dpi), with i = 1,2, 3. In order 
to describe the state of the system in the /i-space, we divide it into celis correspond- 
ing to the macroscopic uncertainties Sqt and Spi, each one occupying in the \i -space a 
volume 



These celis must be sufïiciently small related to the macroscopically measurable 
dimensions but also sufïiciently large to contain a great number of partides. 

When considering the system as a whole, its macroscopic state is given by a collection 
of points rii corresponding to the number of partides inside each celi. Now, if we take 
into account the 6N-dimensional phase space V, in which each of the states assumed by 



Wi = óq! 5q 2 Sq 3 ópi 5p 2 5p 3 . 
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the system is represented by a point, to each configuration m corresponds in T a celi 
with volume given by 



N 



(sv) r = n (wt) 



n=l 



Considering that the permutation of the partides inside the celis of the T space does 
not change the macroscopic state of the system, then to each collection of numbers rii 
corresponds a volume Í2 n in the T space B given by 



The state of the system is determined by the star occupied by the representative 
point of the system in the T space. This way, macroscopically, it is only possible to 
distinguish in which star the system is, such that any point in this star corresponds 
to a same macroscopic state. When we consider a system which is not in equilibrium, 
a change in its macroscopic state can only be observed when the point describing the 
system changes star. The crossing time is small but finite. During this period of time 
the macroscopic state of the system does not change notwithstanding its microscopic 
state is continuously changing. 

Thus, from the point of view of statistical physics, the introduction of a fundamental 
interval of time appears in a very natural way. That is still more significant when we 
remember that the predictions of QM are always obtained as mean vàlues of observ- 
ables. The uncertainty relations, according to the usual intèrpret ation of QM — the 
Copenhagen interpretation — , are independent of the arguments above. If we accept 
that they play a fundamental role in the microscopic world — and this is postulated by 
Copenhagen — , then the concept of chronon, as a fundamental interval of time, must 
be related to them. 

5.2 Discretized Liouville equation and 

the time— energy uncertainty relat ion 

An attempt to set up a relationship between the chronon and the time-energy un- 
certainty relation has been put forward by Bonifacio (1983) |72j, extending the coarse 
graining hypothesis to the time coordinate. In the conventional QM the density oper- 
ator evolves according to the Liouville-von Neumann equation 

.lanecl 71 calls it a star. 
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^ = -iCp(t) = -í [È,p], (99) 

where £ is the Liouville operator. One can immediately observe that, if H is time 
independent, the solution is given by 

p(T) = exp (-i^j p(0) exp (i^j , (100) 

which gives the time evolution of the density operator starting from an initial time to, 
such that T = t — í is the evolution time. 

When we build a coarse grained description of the time evolution, by intro- 
ducing a graining of value r such that the evolution time is now given by T = 
kr (k = 1, 2, . . . , oo), we have that the resulting density operator p does not satisfy 
the continuous equation fl99l) but a discretized form of it given by 

mzBzJl = - icm , (ioi) 

T 

with t = kr, which reduces to the Liouville- von Neumann equation when r — > 0. In 
the energy representation |n), once satisfied certain conditions which ensure that p{k) 
is a density operator, we have that fllOip rules for p an evolution which preserves trace, 
obeys the semigroup law and is an irreversible evolution towards a stationary diagonal 
form. In other words, we observe a reduction of state in the same sense as in the 
measurement problem of QM. This reduction is not instantaneous and depends on the 
characteristic value r: 

p(t) *^°^p nri (0)|n)(n| . 



It is important to observe that the non diagonal terms tend exponentially to zero 
according to a factor which, to the first order, is given by 



-^l m Tt 



exp 

Thus, the reduction to the diagonal form occurs provided we have a finite value for r, no 
matter how small, and provided we do not have uj nm r «C 1 for every n,m, where u nm = 
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(E n — E m )/h are the transition freqüències between the different energy eigenstates. 
This latter condition is always satisfied for systems not bounded. 

These results, together with an analysis of the discrete Heisenberg equation defined 
in terms of the average vàlues of observables 

À(t) = Tr (p(t) A) 

in the coarse grained description, suggest an interpretation of r in terms of the uncer- 
tainty relation AEAt > h/2 such that r is a characteristic interval of time satisfying 
the inequality 

t - Te ^ïXe with AE = ^ {H2) ~ {H)2 ' (102) 

so that the mathematical meaning of the time-energy uncertainty relation is that of fix- 
ing a lower limit for the time interval within which the time evolution can be described. 
Thus, ". . . the coarse grained irreversibility would become a necessary consequence of 
an intrinsic impossibility to give an instantaneous description of time evolution due to 
the time-energy uncertainty relation" . 

Since the density operator, in the energy representat ion, tends to a diagonal form, 
it seems to be tempting to apply it to the measurement problem. We can also observe 
that, even without assuming any coarse graining of time, namely, without using the 
statistical approach adopted by Bonifacio, the reduction to a diagonal form results 
straightforwardly from the discrete Liouville equation and some asymptotic conditions 
regarding the behaviour of the solution, once satisfied[73] the inequality uj nm r <C 1. 
See also [74"] . 

The crucial point, from which derives both the decay of the non-diagonal terms of 
the density operator and the very discrete Liouville equation, is that the time evolution 
operator obtained from the coarse grained description is not a unitary operator. 
This way, the operator 

V(t = hr,t = 0) = ^— T , (103) 

as all the non-unitary operators, does not preserve the probabilities associated with 
each of the energy eigenstates that make up the expansion of the initial state in that 
basis of eigenstates. We must recali that the appearance of non-unitary time evolution 
operators is not associated with the coarse grained approach only, since they also result 
from the discrete Schròdinger equations. 
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5.3 The measurement problem in quantum 
mechanics 

Let us apply the discrete formalism introduced in the previous Subsection to the mea- 
surement problem. Using a quite general formalization, we can describe the measure- 
ment process taking advantage of the properties observed for the evolution of the density 
operator as determined by the discrete Liouville-von Neumann equation0 

When speaking of measurement, we have to keep in mind that, in the process, an 
object O, of which we want to measure a dynamic variable R, and an apparatus A, 
which is used to perform such measurement, are involved. Let us suppose that R is 
the operator associated with the observable R, with an eigenvalue equation given by 
R\r) = r\r) and defines a complete basis of eigenstates. Thus, considered by itself, any 
possible state of the object can be expanded in this basis: 

l^)o = E c -l r )o- (104) 

r 

As regards the apparatus A, we are interested only in its observable A, whose 
eigenvalues a represent the possible vàlues indicated by a pointer. Besides, let its various 
internal quantum numbers be labelled by an index n. These internal quantum numbers 
are useful to specify a complete basis of eigenvectors associated with the apparatus: 

À \a, n) A — a \a, n) A . (105) 

Now, let us suppose that the apparatus is prepared in an initial state given by 
|0, n)^, i.e., in the initial state the value displayed is zero. The interaction between the 
two systems is introduced by means of the time evolution operator and is such that 
there is a correlation between the value of r and the measure a r . In order to deal with 
the measurement process itself we consider a quite general situation. First of all, let us 
consider the following pure state of the system object + apparatus (O + A): 

\V n ) = \V) \0,n) A . (106) 

The evolution of this state, in the continuous description, using the evolution oper- 
ator, is given by 

Ü(t,t o )\y)a\0,n) A = Y, c r\ a r;r,n) = \¥ n ) (107) 

r 

^We follow closely the description exhibited in Ballantine. 83J 
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which is a coherent superposition of macroscopically distinct eigenstates, each one cor- 
responding to a different measure a r . The great problem for the Copenhagen interpre- 
tation results from the fact that it considers the state \^ l n ) as associated with a single 
system: a pure state provides a complete and exhaustive description of an individual 
system. Thus, the coherent superposition above describes a single system so that, at 
the end of the interaction which settles the measurement, the display should not show 
a well-defined output since Í llü7p describes a system which is a superposition of all its 
possible states. 

However, we know from experience that the apparatus always displays a single value 
as the output of the measurement. It is this disagreement between observation and the 
description provided by the formalism, when interpreted according to Copenhagen, 
which results in the necessity of introducing the postulate of the reduction of the vector 
state 

l*n) — * \a ro ;r ,n) 
where r is the value displayed by the apparatus. 

This fact has been considered by many as being a problem for the usual interpre- 
tation of quantum mechanics.|84, 83J The attempts to find a solution, in the context 
of different interpretations, have been numerous, from the Many- Worlds interpretation, 
proposed by Everett and Wheeler.|82j to the measurement theory by Daneri, Loinger 
and Prosperi, E5] in which the reduction of the quantum state is described as a 
process triggered by the appearance of aleatory phases in the state of the apparatus, 
just because of its interaction with the elementary object. The approach introduced 
here is — by contrast — rather simpler. 

As an initial state in the measurement process, let us consider a mixed state for the 
composite system O + A, 

p^Y.Cu K)« , (los) 



where C n is the probability associated with each of the states Such probability 

is, as in the classical physics, an 'ignorance' probability, i.e., it is not intrinsic to the 
system. In the continuous case, when we apply the time evolution operator to that 
density operator we get a final state given by 

p í = Üp i tfY.C n K)(K = (109) 
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P l = c lx c n^C n {|a ri ;ri,n) (a ri ; n, n|} 



(110) 



so that the presence of non-diagonal terms corresponds to a coherent superposition of 
states. In this case, the postulate of the reduction of the quantum state is connected 
with the non-diagonal terms of the density operator. It is postulated that when a 
measurement is carried out on the system, the non-diagonal terms tend instant aneously 
to zero. Since in the continuous case the time evolution of the state results from 
the application of a unitary operator, which preserves the pure state condition p 2 = 
p, it is impossible to obtain the collapse of the pure state from the action of such 
operator. In the diagonal form the density operator describes an incoherent mixture of 
the eigenstates of A, and the indetermination regarding the output of the measurement 
is a sole consequence of our ignorance about the initial state of the system. 

In the discrete case, which has the time evolution operator given by (1 1 3 j) . with 
the interaction between apparatus and object embedded in the Hamiltonian H, the 
situation is quite different. The main cause of such difference lays in the fact that 
the time evolution operator is not unitary. Let us consider the energy representat ion, 
describing the eigenvalue equation of the Hamiltonian as H\n) = E n \j) so that the 
eigenstates \n) are the states with defined energy. From the formalism of the density 
matrices we know that when the operator R is diagonal in the energy representation 
then, when calculating the expected value of the observable, we do not obtain the 
interference terms describing the quantum beats typical of a coherent superposition of 
the states \n). 

As the time evolution operator is a function of the Hamiltonian and, therefore, 
commutes with it, the basis of the energy eigenstates is also a basis for this operator. 
We can now use a procedure identical to the one applied by Bonifacio, and consider 
the evolution of the system in this representation. Thus, we have that the operator 
V (t = kr, t = 0) takes the initial density operator p 1 to a final state for which the non- 
diagonal terms decay exponentially with time: 



(r\V(t = kr,t = 0) \s) 



Pr. 



(111) 



P, 



(1 + %ijJ rs T) 



t/r ' 



with 



1 



(E r - E a ) = -AE,, 



(112) 



LO, 



rs 



67 



The expression (111 ip can be written as 



p rs (t) = p rs (0)e~ Jrst e- lUrst (113) 

such that, 

7rs = ^ln(l + ^ 2 s r 2 ), (114) 

v rs = — tan -1 (u ra r). (115) 

T 

We can observe right away that the non-diagonal terms tend to zero with time and 
the decay is faster the larger the value of r, which here is an interval of time related 
to the whole system O + A. If we keep in mind that in the coarse grained description 
the value of the time interval r originates from the impossibility to distinguish between 
two different states of the system, we must remember that the system O + A is not an 
absolutely quantum system. That means that r could be sígnificantly larger, implying 
an extremely faster dampíng of the non-diagonal terms of the density operator (figure 
15.31) . We then arrived at a process like the one of the reduction of the quantum state, 
even if in a very elementary formalization. This result seems to be very encouraging 
regarding future researches on such important and controversial subject. 

Some points must be pointed out from this brief approach of the measurement 
problem. First of all, we must emphasize that this result does not occur when we use the 
time evolution operators obtained directly from the retarded Schròdinger equation. The 
dissipative character of that equation causes the norm of the state vector to decay with 
time, leading also to a non-unitary evolution operator. However, this operator is such 
that, in the defmition of the density operator we get damping terms which are effective 
even for the diagonal terms. This point, as well as the question of the compatibility 
between Schròdinger's picture and the formalism of the density matrix are going to be 
analysed in Appendix [A] As the composite system O + A is a complex system, it is 
suitably described by the coarse grained description, so that the understanding of the 
relationship between the two pictures is necessary in order to have a deeper insight on 
the processes involved. 

Notwithstanding the simplicity of the approach we could also observe the intrinsic 
relation between measurement process and irreversibility. The time evolution operator 
V meets the properties of a semigroup, so that it does not necessarily possess an inverse: 
and non-invertible operators are related to irreversible processes. In a measurement 
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Figure 11: Vanishing in time of the non-diagonal terms of the density operator for two 
different vàlues of r. For both cases we have used AE = 4 eV. One gets: (a) Slower 
damping for r = 6.26 x 10~ 24 s; (b) faster damping for r = xlCT 19 s. 



process, in which the object is lost just after the detection, we have an irreversible 



Finally, it is worth mentioning that the measurement problem is controversial even 
regarding its mathematical approach. In the simplified formalization introduced above, 
we did not include any consideration beyond those common to the quantum formalism, 
allowing an as clear as possible individualization of the effects of the introduction of a 
fundamental interval of time in the approach to the problem. 

The introduction of a fundamental interval of time in the description of the measure- 
ment problem makes possible a simple but effective formalization of the state reduction 
process. Such behaviour is only observed for the retarded case. When we take into 
account a symmetric version of the Liouville-von Neumann equation the solution is 
given by 



where the diagonal elements do not change with time and the non-diagonal elements 
have an oscillatory behaviour. This means that the symmetric equation is not suitable 
to describe a measurement process, and this is an important distinction between the 
two descriptions. 

It is important to stress that the retarded case of direct discretization of the 
Liouville-von Neumann equation results in the same equation obtained via the coarse 
grained description. This lead us to the consideration of this equation as the bàsic 
equation to describe complex systems, which is always the case when a measurement 
process is involved. 



In this review we tried to get a better insight into the applicability of the various distinct 
formalisms obtained when performing a discretization of the continuous equations. For 
example, what kind of physical description is provided by the retarded, advanced and 
symmetric versions of the Schròdinger equation? This can be achieved by observing the 
typical behaviour of the solutions obtained for each case and, particularly, attempting 
to the derivation of these equations from Feynman's approach. Then we get that the 



process that could very well be described by an operator such as V. 




6 Conclusions 
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advanced equation describes a system that absorbs energy from the environment. We 
can imagine that, in order to evolve from one instant to another, the system must absorb 
energy, and this could justify the fact that, by using Feynman's approach with the usual 
direction of time, we can only obtain the advanced equation. The propagator depends 
only on the Hamiltonian, being independent of the wave function which describes the 
initial state. So, it describes a transfer of energy to the system. 

The retarded equation is obtained by a time reversion, by an inversion of the di- 
rection of the propagator, i.e., by inverting the flux of energy. The damping factor 
characteristic of the retarded solutions refers to a system continuously releasing energy 
into the environment. Thus, both the retarded and the advanced equations describe 
open systems. 

Finally, the symmetric equation describes a system in an energy equilibrium with the 
environment. Thus, the only way to obtain stationary states is by using the symmetric 
equation. 

Regarding the nature of such an energy, it can be related to the very evolution of 
the system. It can be argued that a macroscopic time evolution is possible only if there 
is some energy flux between system and environment. The states described by the sym- 
metric equation are basically equilibrium states, without nett dissipation or absorption 
of energy by the system as a whole. We can also conceive the symmetric equation as 
describing a closed system, which does not exchange energy with the external world. 

On the other hand, when a comparison is made with the classical approach, we can 
speculate that the symmetric equation ceases to be vàlid when the interaction with 
the environment changes fastly within a chronon of time. Thus, phenomena such as 
the collision of highly energètic partides require the application of the advanced or 
retarded equations. The decay of the norm associated with the vector states described 
by the retarded equation would indicate the very decay of the system, i.e., of a system 
abandoning its initial "equilibrium state". The behaviour of the advanced equation 
would indicate the transition of the system to its final state. This speculation suggests 
another interpretation, closer to the quantum spirit. We could consider the possible 
behaviour of the system as being described by all the three equations. However, the 
ordinary QM works with averages over ensembles, which is a description of an ideal, 
purely mathematical reality. The question is that, if we accept the ergodic hypothe- 
sis, such averages over ensembles are equivalent to averages over time. Anyway, the 
quantum formalism always deals with average vàlues when tackling with the real world. 
When the potentials involved vary slowly with respect to the value of the chronon of 
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the system, which means a long interaction time, we have that the contributions due 
to the transient factors coming from the retarded and advanced equations compensate 
each other, and cancel out. Then, in the average, the system behaves according to 
the symmetric equation. On the contrary, when the potentials vary strongly within 
intervals of time of the order of the chronon, we do not have stationary solutions. The 
discrete formalism describes such a situation by makíng recourse, during the interaction, 
to the transient solutions, which will yield the state of the system after the interaction. 
Afterwards, the system will be described again by a symmetric solution. 

The most conservative quantum interpretation would be that of believing that only 
the symmetric equation describes a quantum system. During the interaction process 
the theory does not provide any description of the system, pointing only to the possible 
states of the system after the transient period. The description of the interaction would 
demand one more ingredient: the knowledge of the interaction process (which would 
imply an additional theoretical development, as, for example, the working out of an 
interaction model). 

Besides the question of the physical meaning of the discretized equations, i.e., of 
the type of physical description underlying it, there is also the question of the time 
evolution of the quantum states. The Schròdinger equations describe the evolution of 
a wave function, with which an amplitude of probability is associated. An analogy 
with the electron theory makes us suppose that this wave function does not react 
instantaneously to the external action, but reacts after an interval of time which is 
characteristic of the described system. In discrete QM, the justification of the non- 
instantaneous reaction comes from the fact that the uncertainty principle prevents a 
reaction arbitrarily close to the action application instant. [68] ES] Such uncertainty 
could be related to the very perturbation caused by the Hamiltonian on the state of 
the system, resulting in an uncertainty relation like the Mandelstam and Tamm fTT] 
time-energy correlation. What we meet is a time evolution in which the 'macroscopic' 
state of the system leaps discontinuously from one instant to the other. Therefore, the 
quantum jumps appear not only in the measurement process, but are an intrinsic aspect 
of the time evolution of the quantum system. The difference, in our case, is that the 
jump does not take the system suddenly out of the quantum state it was endowed with, 
but only determines the evolution of that state. 

Another aspect characteristic of the discrete approach is the existence of an upper 
limit for the eigenvalues of the Hamiltonian of a bounded system. In the description 
of a free particle it has been observed the existence of an upper limit for the energy 
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of the eigenfunctions composing the wave packet which describes the particle, but this 
limit does not imply an upper value for the energy of the particle. The existence 
of this limiting value determines the Hamiltonian eigenvalue spectrum within which a 
normalization condition can hold. Once exceeded that value, a transition to the internal 
excited states of the system takes place. This allowed us, e.g., to obtain the muon as 
an excited internal state of the electron. 

It must be noticed the non-línear character of the relation between energy and 
oscillation frequency of a state, and the fact that the theory is intrinsically non-local, 
as can be confirmed by looking at the discretized equations. It must also be stressed 
that the theory described in this paper is non-relativistic. 

Finally, it must be remembered that the symmetric form of the discrete formalism 
reproduces grosso modo the results of the continuous theory. The effects of the intro- 
duction of a fundamental interval of time are evident in the evolution of the quantum 
systems, but they are - -in general — extremely tiny. There have been some attempts 
to find physical situations in which measurable differences between the two formalisms 
can be observed, but till now with little success. [üü] EZl EB] Maybe this could be af- 
forded by exploiting the consequences of the phase shifts caused by the discretization, 
that we saw in Subsections (14. 2 p and (14.31) . Regarding the justifications for introduc- 
ing a fundamental interval of time, let us for instance recali what Bohr[75] replied to 
the famous 1935 paper by Einstein, Podolski and Rosen[86]: "The extent to which 
an unambiguous meaning can be attributed to such an expression as physical reality 
cannot of course be deduced from a priori philosophical conceptions, but . . . must be 
founded on a direct appeal to experiments and measurements" : Considering time as 
continuous may be regarded as a criticizable philosophical position since, at the level 
of experiments and measurements, nature seems to be discrete. 

More important is to recali that, as already mentioned, the new formalism allows 
not only the description of the stationary states, but also a space-time description 
of transient states: The Retarded Formulation yields a natural quantum theory for 
dissipative systems. It is not without meaning that it leads to a simple solution of the 
measurement problem in QM. Since the present review is still in a preliminary form, 
we shall come back to such interesting problems also elsewhere. 
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APPENDICES 



A Evolution Operators in the Schròdinger 
and Liouville— von Neumann Discrete 
Pictures 

When we think of applying the formalism introduced in the previous Sections to the 
measurement problem, the requirement of the existence of a well-defined evolution 
operator comes out. By well-defined we mean, as in the continuous CclSG, cl unitary 
operator satisfying the properties of a group. 

In the continuous case, when the Hamiltonian is independent of time, the time 
evolution operator has the form 



and is a unitary operator which satisfies the condition that H be hermitian. In the 
continuous case, by definition, every observable is represented by a hermitian operator. 
An operator is unitary when its hermitian conjugate is equal to its inverse, such that 



Another important aspect regarding a unitary operator is related to the probability 
conservation. In other words, if the initial state is normalized to one, it will keep 
its norm for all subsequent times. The evolution operator does not change the norm 
of the states it operates on. Thus, we know beforehand that the evolution operators 
associated with the retarded and advanced discretized Schròdinger equations are not 
unitary operators. 

A.l Evolution operators in the Schròdinger 



For the discretized Schròdinger equation the discrete analogue of the time evolution 
operator can easily be obtained. Let us initially consider the symmetric equation that, 
as already remarked, is the closest to the continuous description. After some algebraical 
handling we get the evolution operator as being 




AU = ÀÀ ] = i. 



picture 
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U (t, t ) = exp 



i(t-t ) . _! 
sm 



h 



so that 



\®{x,t)) = Ú(t,t ) \V(x,t )) = exp 



i(í-ío) . _! /ViT 



Thus, if the eigenvalue equation of the Hamiltonian is given by 



H\*(x,to)) = E\*(x,t )) 



we have that 



I^OMo)}- 



|*(x,í)) = exp 



sm Ix) 



I^OMo)) 



As i7 is a hermitian operator, the evolution operator for the symmetric equation 
is also hermitian. However, the existence of a limit for the possible vàlues of the 
eigenvalues of H implies that, beyond such threshold, the evolution operator is not 
hermitian anymore. In fact, if we consider that beyond the threshold the operator H 
has the form 



H — v + ik, 

where u and k are hermitian operators, we obtain in the continuous approach the same 
results obtained in the discrete case. One of the characteristics of a non-hermitian 
operator is just the fact that it does not conserve the norm of the state it acts on. 
For the retarded equation, the evolution operator is given by 



C/(Mo) 

such that, in the limit r — * 0, 



n -(í-Í0)/T 



1 + -TH 

n 



(116) 



lim 



1 + -TH 

n 



-(í-ío)/t 



which is an expression known as the Trotter equality. Taking the conjugate hermitian 
operator W we can verify that this operator is not unitary. In the basis of eigenstates 
of H we can verify that 
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.2 77121 -(t-*Q)/T 



(n|C/ t ?7|n)= 1 + 



n 



is not equal to 1. Thus, that is the reason why the probabilities are not conserved 
for the solutions of the retarded equation. Besides, as the evolution operator for the 
advanced equation is given by 



it can be verified that the formal equivalence between the two equations is obtained 
by the inversion of the time direction and of the sign of the energy. In the relativistic 
case, this is understandable if we remember that, if a transformation changes the sign 
of the time component of a coordinate 4-vector, then it also changes the sign of the 
energy, which is the corresponding element of the energy-momentum 4-vector. Then the 
retarded equation describes a particle endowed with positive energy travelling forward 
in time, and the advanced equation describes an object with negative energy travelling 
backwards in time, i.e., an antiparticle. [79| I8ÜI |8T| [78] 

A.2 Evolution Operator in the density matrix 
picture 

For the sake of simplicity let \ip(t)) be a pure state. The density of states operator is 
defined as 



It can be shown that such operator evolves according to the following dynamic laws. 
For the retarded case, 



Ü(t,t )= 1 




-, (t-ío)/r 



p(t) = \m)(mi 



Arp(í) 



1 TA 



H(t),p(t)]--.H(t)p(t)H(t); 



for the advanced case 



/W(t) = s [■ 

and, finally, for the symmetric case, 



H(t),p(t)]+-^H(t)p(t)H(t); 



Ap(t) 



1 r 



H(t),p(t) 



ih 
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We can thus observe that the retarded and the advanced equations cannot be ob- 
tained by a direct discretization of the continuous Liouville-von Neumann equation. 
Such formal equivalence occurs only for the symmetric case. Taking into account the 
retarded case, we can obtain the equivalent time evolution operator as being 



V(t,t Q ) 



1 



1 + 



H 



1 (t-ío)r • 



(117) 



We must remark that this operator is different from the one obtained from the coarse 
grained approach, 



Vcg(Mo) 



(t-to)r ■ 



(118) 



and that it is not unitary as well. Quantity Vcg is so defined to have the properties 
of a semigroup: without having necessarily an inverse, but possessing the other group 
properties such as commutativity and existence of an identity (besides the translational 
invariance of the initial condition). 

What we can conclude from such a difference between the two operators is that, 
apparently, the descriptions clash with each other. In the coarse grained approach the 
starting point was the continuous Liouville-von Neumann equation and, by introducing 
the graining of the time coordinate, an evolution operator was obtained satisfying the 
retarded equation 



Arp(í) 



1 

ih 



H(t),p(t) 



The second path consisted in starting from the definition of the density operator, in 
order to determine the dynamical equation it satisfies, and then obtaining the evolution 
operator. 

For the symmetric case the evolution operator is given by 



V(t,t 



o 



exp 



i(t - t ) . _ x (tL 



sm 



t V h , 

which is similar to the operator obtained for the continuous case. 



(119) 



A.3 Compatibility between the previous pictures 

We thus have two distinct evolution operators for the retarded Schròdinger and Li- 
ouville equations so that, once established a connection between them, we arrive at 
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the compatibility of the two descriptions. We try to set up a relation between those 
operators by observing their action on the density operator. So, we expect that both 
operators satisfy the expression 



V(t, t - r)po = Ü(t, t - T)po&(t, t-r), 

where the different action of the operators is basically due to the bilinearity of the 
operator V given by fll 17j> . while U, given by ( 11161) . is linear. This relation is vàlid in 
the continuous case, where the evolution operators act on the density operator according 
to 



pit) = exp [— iC(t — t )/h]p = exp [—iH(t — t )/h]po exp [iH(t — t )/h}. 
Considering the basis of hamiltonian eigenstates |n), we have 



(n\£p(0)\m) = (E n - E m )p nm (0), 



so that 



exp (-iCt)p(0) = exp [-it(E n - E m )]p nm (Q), 



(120) 



exp (-iHt)p(0) exp (iHt) = exp [-it(E n - E m )}p nm (0). 



(121) 



The question is to know if the same is vàlid for the discrete case. For the retarded 
approach we must check whether the relation 



1+ 1 ÍC+ T W H...H 



(í-ío)/r 



Po 



-tH 



(t-to)/T 



1 



WH 



(í-í )/t 



is vàlid. What we get is that, if we consider equations such as f i 1 2 j) and (j!2ip to 
continue to be vàlid in the discrete case, then the above relation is vàlid. For a genèric 
element of the operator we then get 



1 + jr(E n — E m ) + ^E n E n 



JJÏ Pnm(0) 



1 



1 + irE n 



í/r P nm 

(0) 



1 - ïrE m 



t/T' 
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Such equivalence can be also observed for the other cases. However, when we con- 
sider the evolution operator obtained from the coarse grained approach we find an 
incompatibility with the operator deriving from the Schròdinger one. For the operator 
(11161) we have 



The question is to know what is the fundamental difference between the two descrip- 
tions: if both are vàlid and under what conditions. Some points must be emphasized. 
First of all, we must remember that the coarse grained description is a semi-classical 
approach which assumes a system with a certain degree of complexity, while the vec- 
tor state description is a fundamentally quantum approach without any imposition, in 
principle, on the number of degree of freedom of the system described. There is a bàsic 
difference even in the way of conceiving the chronon. In the coarse grained approach it 
is understood as a magnitude inwardly connected to the experimental limitations or, for 
an ideal measurement device, to the limitations imposed by the uncertainty relations. 
For the Schròdinger equation, the value of the chronon is taken as a fundamental inter- 
val of time associated with interaction processes among the components of the system, 
and of the system as a whole with some external potential; i.e., it is associated with 
the infernal processes of the system (as it has been conceived for the classical electron). 
In this way, the absence of the mixed term in the evolution operator obtained with 
the semi-classical procedure is comprehensible, as well as its incompatibility with the 
pure quantum description provided by the Schròdinger equations. As a semi-classical 
approach, the range of applicability of the coarse grained formalism extends to the cases 
where the system to be studied is not purely microscopic, particularly in the measure- 
ment processes. We have to stress the fact that, in this formalization, only the retarded 
equation was obtained. Thus, the system described dissipates energy, i.e., it is an open 
system. This is just the characteristic that allows us to have access to the output of a 
measurement. 

In connection with the operator obtained directly in the Schròdinger picture for the 
retarded case, we have that all the elements of the density matrix, even the diagonal 
ones, are damped with time. Besides, there is also the controversy linked to the non- 
existence, in QM, of an applicability limit of the theory due to the number of degrees 
of freedom involved. The formalism does not distinguish between a microscopic and 
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a macroscopic system, so that it should reproduce what is obtained with the coarse 
grained formalism. This means that the measurement problem also appears in the 
discrete formalism through the non-equivalence of the evolution operators (11171) and 



B Non-Hermitian Operators in the Discrete 
Formalism 

One of the features we have been stressing through this work is the non-hermitian 
character of the discrete formalism. In the Schròdinger representation, for example, the 
continuous equation can reproduce the outputs obtained with the discretized equations 
once we replace the conventional Hamiltonian by a suitable non-hermitian Hamiltonian 
we have called 'equivalent Hamiltonian'. One of the characteristics of a non-hermitian 
operator is that its eigenvalues are defined over the complex number field. A linear 
non-hermitian operator can always be considered as being composed by a hermitian 
part, which supplies the real component of the eigenvalues, and by an anti-hermitian 
part, which gives the complex component. 

In the continuous case, let us take the Hamiltonian as being a non-hermitian operator 
given by 

H = v + ik, 

where u and k are hermitian. Then we have, in the Schròdinger picture, that the time 
evolution operator is given by 



For the discrete case, we get from Appendix |A] that the evolution operator for the 
retarded states is given by (11161) 



íusd. 




n 



(122) 




n -(í-ío)/r 



where H is the hermitian operator associated with the conventional hamiltonian. This 
evolution operator can be written as 



U T et(t, to) 
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cxp 



2r 



exp 



(124) 



Comparing (11221) and (11241) we obtain the equivalence of the hamiltonians once v 
and k are given by 



K 



h , /riï 

z> = — tan 

r \ h 



2r V ' 

For the advanced case we obtain the same expressions except by a minus sign for k. 
For the symmetric case, below the critical limit, we have 



~ -i tH\ 



k = 0. 

Above that limit v ceases to be hermitian and, in this case, the evolution operator 
can be written as 



Ü sym (t,t ) = exp 



so that 



VK 



(t-t c 

2t 



exp 



{t - 1 ) 



ln 



T 



27' 



tH 
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r 



tH 
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\ 




with v being now independent of the Hamiltonian and k ceases to be zero. 

The expressions obtained above show the characteristics that v and k must fulfil in 
order that the continuous equation reproduces the outputs of the discretized equations. 
By observing the continuous evolution operator we have that the anti-hermitian part 
of H shows a non-stationary behaviour, resulting in a damping or amplifying term 
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associated with the evolution of the quantum state. Thus, the stationary solutions 
appear only for the symmetric case below the critical limit. In all the other cases the 
transient term always appears. 

In QM, the non-hermitian operators have been used only as mathematical shortcuts, 
as in the case of the Lippmann-Schwinger equation in the scattering theory. It has 
already been observed that the introduction of such operators could make possible the 
description of unstable states, by phenomenologically linking the transient factor to the 
lifetime of the considered states. [65] If in a certain instant to = the system is in one of 
the eigenstates \n) of the Hamiltonian H then, if such state is unstable, the probability 
of the system to be found in the same state in a later instant t is 

P n (t) = |<n|^|n)| =exp(-í/r L ), 
and that allows us to specify a lifetime t l , for the retarded being 

TL = i l·i l^Eiy ( 125 ) 
and for the symmetric case, above the critical energy, 



21n 





rE n 




h 



Such lifetimes are connected with states that, in the discretized formalism, are 
intrinsically unstable. Only the retarded equation seems to be associated with quantum 
states which decay with time. If that is really vàlid, we have an expression which could 
be used for phenomenologically determining the value of the chronon. Finally, we 
can conclude that the time discret ization brings forth a formalism which, even if only 
hermitian Hamiltonians are involved, is equivalent to the introduction of non-hermitian 
operators in the continuous QM. 
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